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ABSTRACT 


This  thesis  provides  a  comparison  of  the  accuracies  of  two  methods  for  fitting  a 
Weibull  distribution  to  a  set  of  aircraft  engines  time-between-failure  data. 

One  method  used  is  the  Maximum  Likelihood  Method  and  assumes  that  these  engine 
failure  times  are  independent.  The  other  method  is  a  Modified  Method  of  Moments  procedure 
and  uses  the  fact  that  if  time  to  failure  T  has  a  Weibull  distribution  with  scale  parameter  X  and 
shape  parameter  p\  then  Tp  has  an  exponential  distribution  with  scale  parameter  A,p.  The  latter 
method  makes  no  assumption  about  independent  failure  times. 

A  comparison  is  made  from  times  that  are  randomly  generated  with  a  program.  The 
program  generates  times  in  a  manner  that  resembles  the  way  in  which  engine  failures  occur  in 
the  real  world  for  an  engine  with  three  subsystems.  These  generated  operating  times  between 
failures  for  the  same  engine  are  not  statistically  independent.  This  comparison  was  extended  to 
real  data. 

Although  the  two  methods  gave  good  fits,  the  Maximum  Likelihood  Method  produced 
a  better  fit  than  the  Modified  Method  of  Moments.  Explanations  for  this  fact  are  analyzed  and 
presented  in  the  conclusions. 


THESIS  DISCLAIMER 

The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may  not 
have  been  exercised  for  all  cases  of  interest.  While  effort  has  been  made,  within  the  time 
available,  to  ensure  that  the  programs  are  free  of  computational  and  logic  errors,  they  cannot 
be  considered  validated.  Any  application  of  these  programs  without  additional  verification  is  at 
the  risk  of  the  user. 
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EXECUTIVE  SUMMARY 

This  thesis  provides  a  comparison  of  the  accuracies  of  two  methods  for  fitting  a 
Weibull  Distribution  to  a  set  of  aircraft  engines  time-between-failure  data. 

One  method  used  is  the  Maximum  Likelihood  Method  and  assumes  that  these  engine 
failure  times  are  independent.  The  other  method  is  a  Modified  Method  of  Moments  procedure 
and  uses  the  fact  that  if  time  to  failure  T  has  a  Weibull  Distribution  with  scale  parameter  X  and 
shape  parameter  P,  then  Tp  has  an  exponential  distribution  with  scale  parameter  7^.  The  latter 
method  makes  no  assumption  about  independent  failure  times. 

In  a  fleet  of  fielded  aircraft  engines,  engines  are  repaired  when  a  fault  is  detected.  In 
addition,  some  critical  components  are  inspected/repaired  after  a  designated  number  of  flight 
hours.  When  faults  are  detected  the  engine  is  repaired  and  the  component  is  replaced  on  the 
engine  unless  the  engine  has  accumulated  a  total  number  of  flight  hours  near  to  the  hard 
inspection  time  of  a  critical  component(s).  If  this  is  the  case  the  critical  components(s)  are  also 
inspected/repaired.  In  nearly  all  cases,  the  engines  are  not  completely  overhauled  to  restore 
them  to  new  condition.  After  repair,  components  in  the  engine  have  mixed  ages  with  respect  to 
flight  hours  with  some  components  new  or  repaired.  Consequently  the  times  between  failures 
on  a  specific  engine  are  not  independent. 

A  comparison  is  made  for  times  that  are  randomly  generated  from  a  program.  The 
program  generates  times  in  a  manner  that  resembles  the  way  in  which  engine  failures  occur  in 
the  real  world  for  an  engine  with  three  subsystems.  These  operating  times  between  failures  for 
the  same  engine  are  generated  in  a  way  that  they  are  not  statistically  independent.  This  is 
accomplished  by  using  a  program  in  S-Plus  which  simulates  failure  times  for  an  engine  with 
three  subsystems.  The  failure  time  distributions  of  each  one  of  these  subsystems  are  Weibull 
with  fixed  parameters  determined  a-priori  and  based  on  real  world  experience.  The  simulation 
gathers  three  failure  times  for  each  engine.  The  study  includes  failures  for  two  different  sets  of 
20  and  100  engines  respectively.  Once  this  comparison  was  completed,  the  study  was  extended 


Xlll 


to  real  data  using  TF-34  engine  operating  times  between  failures  from  the  Naval  Aviation 
Logistics  Data  Analysis  (NALDA)  data  base. 

The  estimation  of  the  parameters  for  the  Weibull  Distribution  which  fits  the  failure 
times  was  done  using  spread-sheet  software  programs.  The  accuracy  of  each  method  was 
measured  using  the  Kolmogorov-Smirnov  test  statistic,  embedded  in  the  S-Plus  software 
package. 

Although  the  two  methods  gave  good  fits,  the  Maximum  Likelihood  Method  produced 
a  better  fit  more  often  than  the  Modified  Method  of  Moments.  Explanations  for  this  fact  are 
analyzed  and  presented  in  the  conclusions. 


xiv 


LIST  OF  SYMBOLS,  ACRONYMS  AND  ABBREVIATIONS 


P  shape  parameter  Weibull  distribution 

A.  reciprocal  of  scale  parameter  Weibull  distribution 

0  scale  parameter  Weibull  distribution 

CDF  Cumulative  Distribution  Function 

ERAP  Engine  Reliability  Analysis  Program 

exp  Exponential  Distribution 

\j\  Jacobian 

K-S  Kolmogorov-Smirnov 

MLE  Maximum  Likelihood  Estimator 

MMM  Modified  Method  of  the  Moments 

n  number  of  failure  times 

NALDA  Naval  Aviation  Logistics  Data  Analysis 

NAWCADPPE    Naval  Air  Warfare  Center  Aircraft  Division  Propulsion  and  Power 
Engineering 


5 

Ti 

u 


statistical  estimate  of  the  standard  deviation 

failure  time 

statistical  estimate  of  the  mean 


xv 


XVI 


I.  INTRODUCTION 

A.         BACKGROUND 

1.  Engine  Reliability  Analysis  Program  (ERAP) 

The  Naval  Air  Warfare  Center  Aircraft  Division  Propulsion  and  Power  Engineering 
(NAWCADPPE)  was  tasked  with  analyzing  aircraft  engine  failure  data  to  estimate  aircraft 
engine  reliability  and  maintainability  factors  and  to  assess  the  effectiveness  of  existing  naval 
aircraft  engine  maintenance  policies.  The  ERAP  program  was  established  to  help  accomplish 
this  task.  It  generates  nine  reports  designed  to  identify  problem  areas  and  trends  in  engine 
logistic  support,  time  between  engine  failures,  engine  availability,  engine  readiness  degraders, 
and  other  engine  performance  factors.  The  data  base  used  in  the  ERAP  program  is  the  Naval 
Aviation  Logistics  Data  Analysis  (NALDA)  data  base.  These  reports  provide  summary 
measures  such  as  mean  flight  hours  between  engine  failure  and  other  averages  and  trend  lines 
on  these  averages.  They  cannot  be  used,  for  example,  to  compute  residual  reliabilities  of 
specified  additional  lifetimes  for  engines  with  known  number  of  repairs  and  operating  times 
since  last  failure.  Nor  can  they  be  used  to  determine  optimal  maintenance/replacement  policies 
for  critical  engine  components.  A  more  detailed  analysis  of  the  engine  failure  data  in  the 
NALDA  data  base  is  required  to  make  these  types  of  assessments.  Michael  R.  Caudill  [Ref  1] 
and  John  A.  Malsbury  [Ref.  2]  made  more  detailed  analyses  of  TF-34  engine  performance  in 
their  Master  of  Science  theses  than  those  provided  under  the  ERAP  program.  Mark  E.  Mlikan 
[Ref.  3]  developed  an  integral  equation  for  determining  optimal  no-build  times  for  critical 
components  of  TF-34  engines  relative  to  their  hard  inspection  times.  All  three  of  these  theses 
used  a  Weibull  distribution  to  model  the  probability  distribution  of  the  flight  hours  between 
engine  failures  for  the  TF-34  aircraft  engines  archived  in  the  NALDA  data  base. 


2.  NALDA  Data  Base 

NALDA  is  an  operational  automated  information  system.  It  is  a  centralized  Naval 
Aviation  3M  data  base  supported  by  the  Defense  Megacenter,  Mechanicsburg,  Pa.  A  user  of 
the  NALDA  data  base  must  make  coded  queries  to  the  system  to  obtain  specific  information. 
This  requires  considerable  experience  with  NALDA.  NALDA  contains  numerous  data 
elements  on  Navy  aircraft  engines  by  engine  type  and  serial  number.  Anytime  an  aircraft 
engine  is  inspected,  maintained,  removed  from  an  aircraft,  repaired,  transported,  reinstalled  or 
discarded,  a  report  is  submitted  that  is  translated  into  the  NALDA  data  base.  The  NALDA 
data  base  must  be  "cleansed"  before  it  can  be  used  to  perform  some  detailed  statistical 
analyses  such  as  fitting  a  probability  distribution  to  the  empirical  cumulative  probability 
distribution  of  operating  times  between  engine  failures.  M.  Caudill  did  this  cleansing  of  the 
NALDA  data  for  TF-34  engines  provided  to  him  by  support  staff  at  NAWCADPPE.  Details 
of  this  cleansing  process  and  its  purpose  are  described  in  CaudilTs  thesis  [Ref  1]. 

3.  Operating  Times  Between  Engine  Failures 

The  NALDA  data  base  reports  the  operating  time  between  engine  removals.  Engines 
can  be  removed  for  non-failure  reasons.  Caudill  [Ref.  1]  established  a  data  base  of  operating 
time  between  engine  failures  for  TF-34  engines  by  engine  serial  number.  These  times  are 
"not"  statistically  independent  because  of  the  manner  in  which  aircraft  engines  are  repaired. 
More  often  than  not,  they  are  not  totally  overhauled  before  being  reinstalled.  Only  minor  repair 
is  made  on  some  occasions.  The  operating  time  between  the  first  and  second  failure  of  an 
engine  is  dependent  on  the  operating  time  to  first  failure  and  the  extent  of  the  repair.  Engines 
with  long  first  failure  times  and  minor  first  repair  will  contain  many  used  components  with  large 
accumulated  operating  hours  when  reinstalled.  They  will  likely  have  shorter  operating  times  to 
the  next  failure  than  engines  with  small  operating  times  to  first  failure.  Consequently,  the 
operating  times  between  engine  failures,  for  the  same  engine,  are  not  statistically  independent. 


B.  FITTING  THE  PROBABILITY  DISTRIBUTION  OF  OPERATING  TIME 

BETWEEN  FAILURES  FOR  TF-34  ENGINES 

1.  Weibull  Distribution 

Aircraft  engines  wear  out.  Consequently,  any  probability  distribution  used  to  model 
their  time  to  failure  should  reflect  wear-out  characteristics.  The  Weibull  distribution  has  these 
properties.  It  was  used  by  Caudill  [Ref.  1],  Malsbury  [Ref.  2],  and  Mikan  [Ref.  3]  to  model 
the  distribution  of  operating  times  between  failures  for  TF-34  engines.  They  used  a  modified 
method  of  moments  procedure  to  compute  the  parameters  of  the  fitted  Weibull  distribution. 
This  method  used  special  properties  of  the  Weibull  distribution  and  the  exponential  distribution; 
the  sample  mean  and  the  sample  standard  deviation  were  used  to  estimate  the  parameters  of  the 
fitted  Weibull  distribution.  It  was  chosen  because  it  was  known  to  be  less  sensitive  to  large 
positive  or  negative  correlations  (dependency)  in  the  data  than  the  more  standard  methods  of 
maximum  likelihood  estimation  when  statistical  independence  is  assumed. 

2.  Research  Question 

The  analysis  performed  in  this  thesis  provides  a  comparison  of  the  accuracies  of  the 
maximum  likelihood  estimation  procedure  (assuming  independence)  and  the  modified  method 
of  moments  procedure  used  by  Caudill  for  fitting  a  Weibull  distribution  to  the  probability 
distribution  of  operating  times  between  engine  failures.  These  comparisons  are  made  for  times 
that  are  randomly  generated  in  a  manner  that  resembles  the  way  in  which  engine  failure  times 
occur  in  the  real  world.  These  randomly  generated  operating  times  between  failures  for  the 
same  engine  will  not  be  statistically  independent.  Once  this  comparison  was  completed,  it  was 
extended  to  real  data  using  TF-34  engine  operating  times  between  failures  from  the  NALDA 
data  base.  [Ref  1] 


n.  GENERATION  OF  FAILURE  TIMES 


A.         PROCEDURE 

To  generate  the  failure  times,  a  program  in  S-Plus  (See  Appendix  A)  was  developed 
with  the  following  properties. 

1 .  Each  engine  has  3  components  (subsystems). 

2.  The  failure  time  distributions  of  each  of  these  components  are  Weibull  distributions 
with  fixed  parameters  determined  a-priori  and  based  on  experience. 

3.  The  data  generated  is  a  Type  II  Censored  Data  Set,  and  a  total  of  3  failure  times 
are  gathered  for  each  engine. 

With  these  criteria,  the  program  generates  the  failure  times  in  the  following  way: 

Step  1:  Three  Weibull  random  numbers  are  generated.  They  represent  the  first  failure 
time  that  each  one  of  the  three  components  could  have.  These  values  constitute  the 
vector  faults  1.  Three  additional  Weibull  random  numbers  are  generated.  They 
represent  the  second  failure  time  that  each  one  of  the  three  components  could  have. 
These  values  constitute  the  vector  faults2.  Three  additional  Weibull  random  numbers 
are  generated.  They  represent  the  third  failure  time  that  each  one  of  the  three 
components  could  have.  These  values  constitute  the  vector  faults3. 

Step  2:  The  minimum  value  of  the  vector  faults 1  is  selected,  this  value  is  the  time  of 
the  first  failure  of  the  engine.  The  program  creates  a  new  vector  ?iewfail2  that  is  the 
difference  between  the  vector  faults  J  and  the  minimum  value  selected;  it  is  clear  that 
this  new  vector  has  a  component  that  is  0.  This  component  is  replaced  in  the  vector 
newfail2  by  the  element  of  the  vector  faults2  that  has  the  same  index  as  the  first  failure. 

Step  3:  At  this  point  the  engine  has  two  components  that  have  accumulated  working 
hours  and  one  that  has  been  just  repaired.  The  vector  newfail2  has  two  elements  that 
have  the  remainder  of  the  theoretical  first  failure  for  these  components  and  one  element 
that  is  the  theoretical  second  failure  for  the  other  component.  So,  again,  the  minimum 
element  of  this  vector  is  selected  and  it  is  the  second  failure  time  of  the  engine. 


Step  4:  The  program  creates  a  new  vector  called  newfaiB  that  is  the  difference 
between  the  vector  newfail2  and  this  minimum  value  selected;  it  is  clear  that  this  vector 
has  a  component  that  is  0. 

If  the  second  engine  failure  time  has  taken  place  in  the  same  component  that  had  the  first 
failure,  the  0  element  is  replaced  by  the  element  of  the  vector  faults3  whose  index  matches  that 
of  the  component  that  has  had  the  two  first  failures  or  by  the  element  of  the  vector  faults! 
whose  index  matches  with  the  component  that  had  the  second  failure.  Again,  there  is  an 
element  in  the  vector  faults3  that  is  a  theoretical  failure  time  for  one  component  just  repaired, 
and  two  elements  that  have  the  remainder  of  theoretical  failures.  So,  once  again  the  minimum 
element  of  this  vector  will  be  the  value  of  the  third  engine  failure  time.  The  vector  engine/ails 
gives  the  three  random  failure  times  of  the  engine. 

B.         GENERATION 

An  S-Plus  program  was  developed  to  generate,  in  the  bank  of  engines,  the  three  failure 
times  for  the  engine.  The  program  can  be  run  for  different  numbers  of  engines.  In  this  thesis  the 
program  was  run  for  100  engines  with  a  total  of  300  failure  times  and  for  20  engines  with  60 
failure  times.  This  was  done  to  determine  if  the  number  of  engines  has  an  apparent  effect  in  the 
results  of  the  study.  For  each  one  of  these  engine  sets,  30  replications  were  generated.  These 
replications  were  exported  from  S-Plus  to  Excel  with  the  export  tool  via  the  S-Plus  file-Export 
menu  and  used  to  assess  averages  on  the  accuracy  of  the  two  methods  of  estimation- 

The  first  replication  for  both  20  and  100  engines  can  be  seen  in  the  Appendix  A,  all  of 
the  other  replications  have  been  retained  on  diskettes  attached  to  the  master  copy  of  this  thesis. 


ffl.  ESTIMATION  OF  PARAMETERS  WITH  THE  MAXIMUM  LIKELIHOOD 
AND  THE  MODIFIED  MOMENT  METHOD 


A.         PROCEDURE 

The  two  parameter  Weibull  distribution  analyzed  in  this  thesis  has  a  probability  density 
function  given  by: 

fT(l)  =  ^'"-,e-w',  (1) 

or  equivalently 

p 
fT(t)  =  pe-pt^e-{^   ,  (2) 

1 
where  0  =  —    Customarily  (3  is  called  the  shape  parameter  and  X  (or  9  )  is  called  the  scale 

parameter. 

To  estimate  the  parameters,  the  capabilities  of  the  spread  sheet  program  Excel  were 
used.  The  failure  times  are  generated  in  S-Plus  and  exported  to  Excel,  where  they  are  displayed 
automatically  on  the  Sheet  1  Worksheet  of  Excel. 

The  Maximum  Likelihood  estimators  are  calculated  based  on  well  known  formulas 

[Ref.  4].  Solve  numerically  for  (3 : 

n 

n^Tip  LnTi 
-  +  tL"Tl-—« =°.  (3) 

P    i-i  y  fjfi 


1=1 


Then: 


x  = 


n 


(4) 


0=1. 


(5) 


The  parameters  for  the  Modified  Method  of  Moments  were  calculated  based  on  the 
fact  that  if  T  has  a  Weibull  distribution  then  Tp  has  a  Exponential  distribution  with  a  scale 
parameter  of  A.p .  The  following  demonstration  is  from  [Ref  5]: 

Considering  the  Weibull  distribution  given  by: 


_  RlP/P-^-(^)P 


f,(t)  =  Wt*-le 


(6) 


Let  x  =  tp    so   ?  =  x1/p, 


Then 


dt       1    g-i 
dx      |3 


(7) 


fx(*)  =  M/('(*))  = 


i^1 
P 


/,(*') 


i     I-i  ^ 

P 


(8) 


Therefore. 


X~exp(Zfl) 


(9) 


It  is  well  known  that  the  mean  and  standard  deviation  of  the  exponential  distribution 
are  equal.  Therefore  their  statistical  estimates  u  and  s  should  be  nearly  equal.  That  is 


s  = 


1 


n-\ 


Z(^-(^)): 


/=! 


(10) 


where  the  Tt  are  the  times  between  engine  failures. 


The  solution  for  P  in  equation  (10)  will  be  denoted  by  (3 

p 


Since  E(T,P)  =  Q 


P\_aP 


IV  1   " 

,  the  sample  mean,  —  ]T  J,73  estimates  9  p ,  therefore. 


0  = 


1    " 
Lw£? 


(11) 


This  is  the  same  estimate  obtained  for  the  MLE  of  0. 


B. 


ESTIMATION 


To  carry  out  the  estimation  of  parameters  with  the  two  methods,  two  Macros  in  Visual 
Basic  were  developed.  One  Macro  for  20  engines  (See  Appendix  B),  that  uses  60  failure  times, 
and  the  other  for  100  engines  (See  Appendix  C),  that  uses  300  failure  times. 

Basically,  these  two  Macros  are  similar  and  therefore  they,  and  procedures  for  using 
them,  can  be  similarly  described.  Once  the  data  is  exported  to  Excel,  it  appears  on  Sheet  J.  At 
this  point  it  is  necessary  to  open  the  file  with  the  Macro,  called  Macro20.xls  for  the  sample  of 
20  engines  and  Macrol00.xls  for  the  sample  of  100  engines  respectively.  Next  the  window 
with  the  Macro  is  hidden  from  the  same  Sheet  /.Then  the  Macro  can  be  run.  To  run  the  Macro, 


choose  Macro  from  the  Tools  menu.  From  Macro/Name  Reference  choose  the  name  of  the 
Macro  and  click  on  run.  The  Macro  does  all  the  previous  calculations  and  operations  needed 
to  estimate  the  parameters  and  to  entitle  the  entire  Worksheet.  After  these  calculations  the 
Macro  stops.  At  this  point  it  is  necessary  to  introduce  the  following  data: 

In  Cell  B6  introduce  the  value  previously  calculated  by  the  macro  for  the  initial 
value  of  the  shape  parameter  which  is  in  Cell  B4. 

This  initial  parameter  was  introduced  to  assure  the  convergence  with  the  Tool  Goal 
Seeker,  because  without  it,  a  solution  may  not  be  reached.  The  initial  value  is  given  by  [Ref. 

6]: 

i 

n  n 

YsiLnTif-fcLnTifln 


Po  = 


6^ 

7f 


i=l  (=1 


n-\ 


2 


(12) 


Using  these  choices  for  the  initial  shape  parameter,  an  average  of  only  3.5  Newton 

iterations  were  needed  to  achieve  four-place  accuracy  in  the  solution  for  p  . 

After  introducing  the  initial  value  given  in  Cell  B4  into  Cell  B6,  cell  C6  must  be 
highlighted.  Then  the  Tools  Menu  is  selected,  Goal  Seek  is  selected,  a  0  is  placed  in  the  "To 
Value,"  and  B6  is  entered  as  the  "Changing  Cell." 

Equation  (3)  is  in  Cell  C6  and  the  Goal  Seeker  searches  for  the  value  of  (3  which 
satisfies  that  equation.  There  were  no  non-convergence  problems  in  any  of  the  70 
replications  made  with  that  formula.  Once  the  solution  is  achieved,  the  final  value  for  the 
shape  parameter  appears  in  Cell  B6.  The  solutions  for  X  and  9  appear  in  Cell  B9 
(Equation  (4))  and  Bl  1  (Equation  (5))  respectively.  These  solutions  for  /?,  X  and  0  are 
the  maximum  likelihood  estimates  for  these  parameters  of  the  Weibull  distribution  that  is 
used  to  model  the  probability  distribution  of  the  family  of  times  between  failures  on  all 
engines.  These  values  are  printed  on  the  left  side  of  Sheet2. 

On  the  right  side  of  Sheet  2  are  the  calculations  carried  out  by  the  Modified 
Moment  Method. 
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There  were  convergence  problems  in  calculating  P  using  Equation  (10).  Amazingly 
the  same  formula  may  converge  or  not  depending  on  the  terms  on  each  side.  Through  trial  and 
error  it  was  discovered  that  by  writing  the  formula  in  the  following  way  it  was  possible  to  get 
convergence  in  all  replications: 


n       _ 


£77' 


V^T 


£(77* -(CT>))2 


L'=i 


3T  (13) 


In  Cell  Q6  it  is  necessary  to  introduce  the  initial  shape  parameter  given  in  Cell  Q4 
(Equation  (12)).  Once  done,  highlight  Cell  P7  (Equation  (13)),  select  the  Tools  Menu,  select 
Goal  Seek  and  depending  on  the  number  of  engines,  enter  in  "To  Value"  7.811334  (20 
engines,  60  failure  times,  60/^59  )  or  17.349447  (100  engines,  300  failure  times,  300/ V299  ). 
In  the  "Changing  Cell,  Q6,"  the  Goal  Seeker  looks  for  the  P  value  that  will  satisfy  Equation 
(13).  Cells  Q9  and  Qll  contain  the  values  for  the  scale  parameter  X  and  6  determined  by 
Equation  (11)  and  (5). 

A  total  of  30  replications  for  each  set  of  engines  were  done.  Afterwards,  for  the 
reasons  explained  later,  10  more  replications  were  done  using  real  data. 

Printouts  of  the  failure  times  and  parameter  calculations  for  the  20  and  100  engines  first 
replication  are  provided  in  Appendix  E. 

A  summary  of  the  parameters  obtained  is  shown  in  the  following  tables: 
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Table  1.  Twenty  Engines  =  60  Failure  Times 

20  Engines  =  60  Failure  Times 

Maximum  Likelihood  Estimator  Modified  Moment  Method 

Shape  Scale  Shape  Scale 


1 

0.973000000 

251 .705300000 

2 

1.105104390 

279  481900000 

3 

1 .052707000 

242.343600000 

4 

1 .005354000 

282.577600000 

5 

1.211857000 

272.413900000 

e 

1.523761000 

277.443800000 

7 

1.053815000 

243.377800000 

$ 

1.015065000 

330.527000000 

9 

0.961137000 

251 .822700000 

1P 

1.108272000 

263.760200000 

11 

1.258861000 

317.378200000 

•- 

12 

1.358459000 

245.953800000 

1? 

1 .323305000 

281 .038800000 

E 

3 

14 

1 .054978965 

223.276800000 

Z 

15 

1.228101693 

264.364500000 

_© 

1? 

1.161814103 

243.734600000 

17 

1.355617000 

268.067900000 

"El 

18 

1 .308393000 

309  037200000 

OS 

19 

1.391697000 

284.440700000 

20 

1.141031467 

242.289400000 

21 

1 .336786000 

271  132800000 

22 

1 .389460000 

244.898600000 

2? 

1 .200747000 

288.213200000 

24 

1 .297427000 

293.585500000 

25 

1.400818000 

310.907900000 

26 

1.231510000 

277.128900000 

27 

1 .048038281 

277.331100000 

28 

1 .093433000 

273.733700000 

,5$ 

1.295133208 

246.789500000 

30 

1  224688831 

303.055100000 

0.985725000 

253.139700000 

1.129925000 

281 .796500000 

1 .036553000 

240.792500000 

1 .079900000 

291 .030800000 

1.102670000 

263.279100000 

1.442500000 

273.090800000 

1 .309674000 

262.946000000 

1.006751000 

329.425200000 

1.078819000 

264.611000000 

1 .084296000 

261 .569300000 

1.395249000 

328.224900000 

1.211490000 

237.176300000 

1.461631834 

289.757400000 

1 .24997981 1 

236.965600000 

1.182801000 

260.845800000 

1 .239576000 

249.537300000 

1.367488000 

268.798900000 

1 .357252000 

312.609300000 

1 .406588000 

285.359900000 

1.324936416 

255.330700000 

1 .223524000 

263.236700000 

1.280865198 

238.618200000 

1.263844000 

293.349300000 

1.260780021 

290.891900000 

1.325062000 

305.455500000 

1.320218000 

283.439100000 

1.043436000 

276.830600000 

1.248177000 

286.965800000 

1.240951783 

243.300600000 

1 .402646081 

317.067700000 

Average 
Std.  Dev 


1 .203679065 
0.145886716 


272.060400000 
25.376732459 


1 .235443638 
0.133749260 


274.848080000 
25.937151577 
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Table  2.  One-Hundred  Engines  =  300  Failure  Times 


100  Engines  =  300  Failure  Times 

Maximum  Likelihood  Estimator  Modified  Moment  Method 

Shape  Scale  Shape  Scale 


1 

1.261028000 

292.556900000 

2 

1.147177000 

258.382300000 

3 

1.192249000 

267.935700000 

4 

1  191393000 

301 .254500000 

5 

1.100167000 

265.867700000 

6 

1.132621000 

273.303500000 

7 

1 .220806000 

282.032200000 

8 

1.206135000 

267.017000000 

9 

1 .254734000 

291 .692500000 

10 

1 .284450000 

285.463200000 

11 

1 .277785000 

285.850800000 

12 

1 .274025000 

263.121900000 

13 

1.192249000 

267.935700000 

u 

14 

1.194163000 

271 .996400000 

at 

— 

15 

1.216995000 

272.494100000 

E 
= 

16 

1.169718000 

293.370600000 

Z. 

M 

1 .283386000 

255.640200000 

_£ 

1? 

1.416971000 

294.799600000 

19 

1  264662000 

268.682800000 

a 

20 

1 .342221 000 

279.712400000 

0! 

21 

1 .229806000 

287.109500000 

22 

1.199120000 

267.635600000 

23 

1.189326000 

288.996500000 

24 

1.199541000 

263.550100000 

25 

1 .268509000 

261.778100000 

26 

1.257310000 

256.670600000 

27 

1 .274659000 

270.850700000 

29 

1 .208075000 

276.652800000 

29 

1 .220756000 

264.01 3900000 

30 

1.127942000 

279.203900000 

1 .288232000 

294.659500000 

1.182793000 

261 .246400000 

1 .250620000 

272.423600000 

1 .243339000 

305.757700000 

1.187671000 

272564100000 

1.205131000 

279.422100000 

1 .334481 000 

290504100000 

1 .246956000 

270.084700000 

1.271828000 

293  041800000 

1.292318000 

286.039200000 

1 .387920000 

293508300000 

1.341512000 

267591700000 

1 .250620000 

272.423600000 

1 .228539000 

274.677200000 

1.231497000 

273.612100000 

1 .200586000 

296.122500000 

1 .259945000 

254.068200000 

1 .422087000 

295.112700000 

1 .379585000 

276.306900000 

1 .448534000 

286  372700000 

1 .269967000 

290.254800000 

1.281756000 

273.804000000 

1.207016000 

290.494900000 

1  346121000 

273.747200000 

1.201983000 

256.990000000 

1 .249969000 

256.167900000 

1.311836000 

273.423600000 

1 .237968000 

278.973700000 

1.219361000 

263.912300000 

1.173691000 

283.295300000 

Average 
Std.  Dev. 


1 .226599300 
0.063780678 


275.185723333 
1 2.379509005 


1.271795400 
0.070516179 


278553426667 
1 3.009908323 
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IV.    ACCURACY  ANALYSIS 

A.  PROCEDURE 

The  accuracy  of  each  method  was  measured  using  the  Kolmogorov-Smirnov  goodness 
of  fit  Test  and  the  S-Plus  software  package.  The  null  hypothesis  in  S-Plus  is  "the  True  cdf 
(empirical)  equals  the  Weibull  distribution".  This  test  has  the  advantage  that  it  uses  continuous 
populations.  Therefore  it  is  not  necessary  to  group  the  data. 

This  test  has  the  disadvantage  that  when  distribution  parameters  are  estimated  from  the 
same  sample  of  data,  the  test  produces  a  conservative  result;  that  is,  the  probability  of  a  Type  I 
error  will  be  smaller  than  specified  [Ref  6].  This  drawback  has  only  a  relative  importance  in  the 
present  study.  The  K-S  Test  p-value  is  a  conservative  result  whether  the  parameters  are 
estimated  by  the  Maximum  Likelihood  Method  or  the  Modified  Moment  Method. 

A  S-Plus  function  called  Tocompare  (See  Appendix  E)  was  created  to  calculate  the  K- 
S  statistics  for  the  set  of  parameters  calculated  in  each  replication  using  both  methods.  This  S- 
Plus  function  computes  the  K-S  Test  Statistic  (Supremum)  and  the  p-value  for  each  method, 
and  it  also  plots  the  empirical  and  hypothesized  distributions.  In  Appendix  E  the  reader  will  find 
the  output  of  the  function  for  the  first  replication  of  each  set  of  parameters. 

B.  ACCURACY 

The  calculations  for  all  replications  and  the  values  obtained  for  each  set  of  engines  are 
presented  in  Tables  3  and  4.  The  letter  "  B"  marks  the  method  that  achieved  a  better  fit,  using 
as  reference  the  K-S  statistic;  that  is,  the  method  with  the  smaller  K-S  statistic  is  the  more 
accurate  procedure  for  that  replication. 

For  the  set  of  20  engines,  the  MLE  method  gave  a  better  fit  in  17  of  the  30  replications. 
For  the  set  of  100  engines,  the  MLE  method  gave  a  better  fit  in  20  of  the  30  replications. 
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Table  3.  Better  Fit  Twenty  Engines 


shape 


20  Engines  =  60  Failures  times 

M.L.E.  M.M.M. 

scale  k-s  p-value  shape  scale  k-s 


p-value 


1 

0.973000000 

251.705300000 

0.054600000 

0.994000000 

0.985725000 

253.139700000 

0.050900000 

0.997700000  B 

2 

1.105104390 

279.481900000 

0.088300000 

0.737300000 

1.129925000 

281.796500000 

0.081000000 

0.826600000  B 

3 

1.052707000 

242.343600000 

0.076300000 

0.876000000 

1.036553000 

240.792500000 

0.073000000 

0.906600000  B 

4 

1 .005354000 

282.577600000 

0.061700000 

0.976300000 

B 

1 .079900000 

291.030800000 

0.086100000 

0.765600000 

5 

1.211857000 

272.413900000 

0.101800000 

0.562800000 

1.102670000 

263.279100000 

0.097300000 

0.621600000  B 

6 

1.523761000 

277  443800000 

0.086700000 

0.757300000 

1 .442500000 

273.090800000 

0.081000000 

0.826400000  B 

7 

1.053815000 

243.377800000 

0.091500000 

0.696500000 

B 

1 .309674000 

262.946000000 

0.115700000 

0.398400000 

6 

1.015065000 

330.527000000 

0.052700000 

0.996200000 

B 

1.006751000 

329.425200000 

0.053100000 

0.995900000 

9 

0.961137000 

251.822700000 

0.088000000 

0.742000000 

B 

1.078819000 

264.611000000 

0.098600000 

0.604700000 

10 

1.108272000 

263.760200000 

0.064500000 

0.964400000 

1 .084296000 

261.569300000 

0.061500000 

0.977100000  B 

11 

1.258861000 

317.378200000 

0.101500000 

0.566300000 

B 

1 .395249000 

328.224900000 

0.126600000 

0.290900000 

12 

1 .358459000 

245.953800000 

0.104300000 

0.531000000 

B 

1.211490000 

237.176300000 

0.121500000 

0.338000000 

13 

1 .323305000 

281 .038800000 

0.075700000 

0.881800000 

B 

1.461631834 

289.757400000 

0.110800000 

0.452600000 

14 

1 .054978965 

223.276800000 

0.088000000 

0.741900000 

1.249979811 

236.965600000 

0.087200000 

0.751700000  B 

15 

1.228101693 

264.364500000 

0.070300000 

0.928000000 

1.182801000 

260.845800000 

0.065700000 

0.958200000  B 

16 

1.161814103 

243.734600000 

0.078800000 

0.850400000 

B 

1 .239576000 

249.537300000 

0.101500000 

0.567200000 

1? 

1.355617000 

268.067900000 

0.070200000 

0.929200000 

B 

1 .367488000 

268.798900000 

0.072500000 

0.910600000 

18 

1 .308393000 

309.037200000 

0.049400000 

0.998600000 

B 

1 .357252000 

312.609300000 

0.054400000 

0.994300000 

19 

1.391697000 

284.440700000 

0.083700000 

0.794200000 

1 .406588000 

285.359900000 

0.083400000 

0.798500000  B 

20 

1.141031467 

242.289400000 

0.091800000 

0.692600000 

B 

1.324936416 

255.330700000 

0.116800000 

0.386000000 

21 

1 .336786000 

271.132800000 

0.107000000 

0.498400000 

1.223524000 

263.236700000 

0.077500000 

0.863800000  B 

22 

1 .389460000 

244.898600000 

0.085000000 

0.778500000 

B 

1.280865198 

238.618200000 

0.092700000 

0.681600000 

23 

1 .200747000 

288.213200000 

0.083700000 

0.795100000 

B 

1 .263844000 

293.349300000 

0.096600000 

0.629800000 

24 

1.297427000 

293.585500000 

0.058400000 

0.986700000 

B 

1.260780021 

290.891900000 

0.059800000 

0.982700000 

25 

1.400818000 

310.907900000 

0.071100000 

0.921800000 

1 .325062000 

305.455500000 

0.071100000 

0.922100000  B 

26 

1.231510000 

277.128900000 

0.062500000 

0.973200000 

1.320218000 

283.439100000 

0.057400000 

0.989000000  B 

27 

1 .048038281 

277.331100000 

0.078500000 

0.853300000 

1.043436000 

276.830600000 

0.077100000 

0.867800000  B 

28 

1 .093433000 

273.733700000 

0.078000000 

0.859200000 

B 

1.248177000 

286.965800000 

0.084100000 

0.790300000 

29 

1.295133208 

246.789500000 

0.055600000 

0.992500000 

B 

1.240951783 

243.300600000 

0.067600000 

0.947000000 

30 

1.224688831 

303.055100000 

0083900000 

0.792200000 

B 

1.402646081 

317.067700000 

0.095900000 

0.639300000 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

1 .203679065 

272.060400000 

0.078116667 

0.822256667 

17 

1.235443638 

274.848080000 

0.083946667 

0.756066667  13 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

0 145886716 

25.376732459 

0.015610639 

0  146004782 

0.133749260 

25.937151577 

0020788358 

0214131410 
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Table  4.  Better  Fit  One-Hundred  Engines 


shape 


100  Engines  =  300  Failures  times 

M.L.E.  M.M.M. 

scale  k-s  p-value  shape  scale  k-s 


p-value 


1 

1.261028000 

292.556900000 

0.031900000 

0.920600000 

B 

1.288232000 

294.659500000 

0.038200000 

0.774000000 

2 

1.147177000 

258.382300000 

0.027400000 

0.977600000 

B 

1.182793000 

261.246400000 

0.032700000 

0.905100000 

3 

1.192249000 

267.935700000 

0.047000000 

0.522300000 

B 

1.250620000 

272.423600000 

0.051000000 

0  416900000 

4 

1.191393000 

301.254500000 

0.055800000 

0.306900000 

1 .243339000 

305.757700000 

0.046600000 

0.532800000  B 

5 

1.100167000 

265.867700000 

0.060600000 

0.220000000 

1.187671000 

272.564100000 

0.045600000 

0.561900000  B 

6 

1.132621000 

273.303500000 

0.043200000 

0.630900000 

1.205131000 

279.422100000 

0.042500000 

0.651400000  B 

7 

1.220806000 

282.032200000 

0.034300000 

0.873100000 

B 

1.334481000 

290.504100000 

0.048900000 

0.469300000 

6 

1.206135000 

267.017000000 

0.037900000 

0.782600000 

1 .246956000 

270.084700000 

0.036300000 

0.823000000  B 

9 

1.254734000 

291 .692500000 

0.021600000 

0.999000000 

B 

1.271828000 

293.041800000 

0.025500000 

0.989900000 

10 

1 .284450000 

285.463200000 

0.022700000 

0.997900000 

1.292318000 

286.039200000 

0.020800000 

0.999500000  B 

11 

1 .277785000 

285.850800000 

0.044600000 

0.589000000 

B 

1 .387920000 

293.508300000 

0.058300000 

0.260200000 

12 

1.274025000 

263.121900000 

0.047400000 

0.511700000 

B 

1.341512000 

267.591700000 

0.048900000 

0.469900000 

13 

1.192249000 

267.935700000 

0.029400000 

0.957700000 

B 

1 .250620000 

272.423600000 

0.035300000 

0.848500000 

14 

1.194163000 

271.996400000 

0.039600000 

0.735300000 

B 

1.228539000 

274.677200000 

0.048600000 

0.477100000 

15 

1.216995000 

272.494100000 

0.031400000 

0.928900000 

1.231497000 

273.612100000 

0.029400000 

0.957500000  B 

16 

1.169718000 

293.370600000 

0.034800000 

0.861700000 

B 

1 .200586000 

296.122500000 

0.041300000 

0.685300000 

17 

1 .283386000 

255.640200000 

0.031800000 

0.922000000 

B 

1 .259945000 

254.068200000 

0.033000000 

0.898900000 

1$ 

1.416971000 

294.799600000 

0.039800000 

0.729300000 

1 .422087000 

295.112700000 

0.038700000 

0.759500000  B 

19 

1 .264662000 

268.682800000 

0.053800000 

0.349300000 

B 

1 .379585000 

276.306900000 

0.065700000 

0.150500000 

20 

1.342221000 

279.712400000 

0.039700000 

0.730800000 

B 

1 .448534000 

286.372700000 

0.067500000 

0.129900000 

21 

1.229806000 

287.109500000 

0.041300000 

0.685900000 

1.269967000 

290.254800000 

0.040300000 

0.713200000  B 

22 

1.199120000 

267.635600000 

0.036000000 

0.830800000 

B 

1.281756000 

273.804000000 

0.044400000 

0.596300000 

23 

1.189326000 

288.996500000 

0.032800000 

0.903200000 

B 

1.207016000 

290.494900000 

0.033600000 

0.886800000 

24 

1.199541000 

263.550100000 

0.032000000 

0.919200000 

B 

1.346121000 

273.747200000 

0.048900000 

0.469200000 

25 

1 .268509000 

261.778100000 

0.028200000 

0.971000000 

B 

1.201983000 

256.990000000 

0.041500000 

0.680900000 

26 

1.257310000 

256.670600000 

0.026400000 

0.984800000 

B 

1 .249969000 

256.167900000 

0.027700000 

0.975800000 

27 

1.274659000 

270.850700000 

0.038700000 

0.759500000 

B 

1.311836000 

273.423600000 

0.048400000 

0482200000 

28 

1 .208075000 

276.652800000 

0.035200000 

0.851300000 

1 .237968000 

278.973700000 

0.028200000 

0.971300000  B 

29 

1.220756000 

264.013900000 

0.033800000 

0.882800000 

B 

1.219361000 

263.912300000 

0.034200000 

0.874100000 

30 

1.127942000 

279.203900000 

0.035700000 

0.838800000 

1.173691000 

283.295300000 

0  031400000 

0.928300000  B 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

1 .226599300 

275.185723333 

0.037160000 

0.772463333 

20 

1.271795400 

278.553426667 

0.041113333 

0.677973333  10 

standev 

standev 

standev 

standev 

standev 

standev 

standev 

standev 

0063780678 

12379509005 

0009086679 

0.209187169 

0.070516179 

13.009908323 

0  010947594 

0.247541499 
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If  the  MLE  and  MMM  methods  are  equally  good  then  one  would  expect  that  the 
probability  the  MLE  gives  a  better  fit  (  using  the  K-S  statistic)  should  be  0.5.  Let  n  represent 
the  probability  the  MLE  gives  a  better  fit;  the  above  sets  of  data  can  then  be  used  to  test  Ho  :  n 
=  0.5  versus  Ha :  n  ^0.5.  For  the  set  of  20  engines,  17  successes  were  observed  out  of  30 
replications  (  MLE  gave  a  smaller  K-S  value  than  MMM);  the  p-value  for  testing  FL :  n  =  0.5 
is  0.5847  so  Ho  would  be  accepted  for  any  a  <  0.5847.  For  the  set  of  100  engines,  20  of  the 
30  replications  gave  successes;  in  this  case  the  p-value  is  0.0987  and  Flo  is  rejected  for  any  a 
>  0.0987.  Certainly  this  comparison  of  MLE  and  MMM,  using  the  binomial  test,  shows  no  real 
difference  for  20  engines  and  marginally  indicates  superiority  of  MLE  for  100  engines. 

If  the  K-S  statistic  mean  and  standard  deviation  are  analyzed  it  is  clear  that  the  MLE. 
gives  the  smaller  mean  and  standard  deviation  of  the  two  cases. 

For  this  data,  it  seems  that  despite  the  fact  that  there  is  some  dependence  between  the 
failure  times,  the  MMM  does  not  give  a  better  fit  than  the  MLE.  This  is  true  at  least  for  the 
failure  times  generated  with  the  simulation  program  used  in  this  study.  This  is  likely  due  to  the 
fact  that  the  failure  times  for  different  engines  (among  engines)  are  independent. 

At  this  point  real  failure  data  was  examined.  This  real  data  is  from  [Ref  1],  and  from 
the  cleansed  data  collected  in  Appendix  B  (Five- Year  Time  Between  Failures  Data  Base).  This 
data  was  cleansed,  removing  zeros  as  failure  times.  Ten  samples  of  data  were  taken  from  the 
data  set  of  operating  times  between  failures  on  TF  34  aircraft  engines.  Five  samples  were  of  60 
failures  times  and  five  samples  were  of  300  failure  times. 

The  same  procedure  to  estimate  the  parameters  and  analyze  the  accuracy  of  the  results 
were  applied.  The  results  are  shown  in  Tables  5  and  6. 

With  the  real  data  the  result  is  the  same,  a  better  fit  is  obtained  more  often  with  the 
MLE  method  than  with  the  MMM  method.  In  the  60  failure  times  case  the  MLE  method  was 
the  better  procedure  in  four  of  the  five  samples,  and  in  the  300  failure  times  case  it  was  better 
three  of  the  five  samples  using  the  K-S  statistic  as  measure  of  accuracy. 
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It  is  important  to  point  out  that  in  these  real  data,  the  engine  structure  is  unknown. 
There  is  no  information  about  how  many  elements  the  engine  has,  or  of  course  the  way  that 
these  components  failed. 

Table  5.  Estimate  Parameters  and  Analyze  Accuracy  Twenty  Engines 


shape 


20  Engines  =  60  Failures  times 

M.L.E.  M.M.M. 

scale  k-s         p-value  shape  scale  k-s         p-value 


1 

1.554927  651.243400  0.103000  0.547400  B 

1.732451   669.495500   0.103200  0.544600 

licati 
Limb 

2 

1.435165  710.748800  0.079700  0.840000  B 

1.735014  748.435600   0.130200  0.261300 

3 

1.531311   603.432400  0.072600  0.909800 

1.635563  613.677100   0.051800  0.997100  B 

gj  ^ 

4 

1.436419  609.475100  0.103800  0.537100  B 

1.485100  615.440400   0.114900  0.406200 

Pti 

5 

1.535733  579.331200  0.135400  0.221300  B 

1  490634  574.988000   0.145900  0.155400 

mean           mean          mean        mean 

mean           mean          mean        mean 

1.414134  642.431620  0.048360  0.511380    4 

1.587754  663.677300   0.057660  0.289300    1 

stdev          stdev          stdev        stdev 

stdev          stdev          stdev        stdev 

0.045711      14  890226  0.008406  0.185231 

0.023623      16.834529  0.006280  0  138620 

C.         ANALYSIS  OF  THE  RESULTS 

The  author  expected  "  a-priori"  a  better  fit  with  the  MMM  than  with  the  MLE,  based 
on  the  fact  that  the  generated  operating  times  between  failures  for  the  same  engine  are  not 
statistically  independent.  But  the  resulting  fit  for  the  data  was  better  with  the  MLE  more  often 
than  with  the  MMM.  This  result  is  further  analyzed  here,  by  examining  how  the  dependence 
between  failure  times  affect  each  method. 

To  assess  dependence  between  failure  times,  the  correlation  and  the  grade  of 
dependence  between  the  failure  times  was  calculated  for  each  replication.  The  correlation 
coefficients  were  calculated  and  added  at  the  end  of  each  replication  (See  Appendix  A).  In 
addition  the  Kendall-Tau  independence  test  was  performed.  The  null  hypothesis  in  S-Plus  for 
this  test  is  "  Tau  is  equal  to  0"  that  is  the  two  samples  are  independent.  This  analysis  was 
performed  for  all  possible  pairwise  combinations  of  the  failure  times  of  an  engine;  namely  first- 
second,  first-third,  second-third  correlations,  tau  coefficients  and  p-values  were  calculated.  For 
the   Kendall-Tau   test   the   minimum,    maximum    and   mean    p-values   of  those   pairwise 
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combinations  of  each  replication  were  added  to  the  tables.  The  results  are  shown  in  Tables  7 
through  10. 

Table  6.  Estimate  Parameters  and  Analyze  Accuracy  One-Hundred  Engines 


shape 


100  Engines  =  300  Failures  times 

M.L.E.  M.M.M. 

scale  k-s         p-value  shape  scale  k-s 


p-value 


.2   <u 

w  S 

~  s 

4>  ^-> 


1 

1.462681   618.916600  0.047200  0.515100 

1.579829  632  033000   0.045700  0.559000    B 

2 

1.439461    663.159800  0.045200   0.571700    B 

1.586168  681.405000   0.063000  0.184300 

3 

1.366409  641.062700  0.040500   0.709200    B 

1.564422  666.063600   0.057300  0.279000 

4 

1.351920  636.911100  0.064600   0.163100 

1.575530  665.621100   0.061900  0.200100    B 

5 

1.450201    652.107900  0.044300   0.597800    B 

1.632821   673.263800   0.060400  0.224100 

mean           mean           mean         mean 

mean           mean           mean         mean 

1.414134  642.431620   0.048360  0.511380    3 

1.587754  663.677300   0.057660  0.289300    2 

stdev          stdev          stdev        stdev 

stdev          stdev          stdev        stdev 

0  045711      14  890226   0.008406   0.185231 

0.023623     16.834529   0  006280  0  138620 

From  those  tables  may  be  deduced  that  with  the  set  of  20  engines  the  correlation  is 
clearly  larger  than  with  100  engines.  When  the  independence  is  tested  with  the  Kendall-Tau 
test,  the  values  of  the  tau's  coefficients  are  similar  to  the  correlation  coefficients,  and  have 
larger  values  for  the  set  of  20  engines  than  the  set  of  100  engines.  In  general  for  this  set  of  data, 
the  correlation  coefficients  are  small  and  the  p-values  in  the  Kendall-Tau  are  not  significant 
enough  to  reject  the  null  hypothesis  of  independence  at  the  usual  values.  That  means  the 
dependence  obtained  in  the  simulation  is  small,  specially  with  the  set  of  1 00  engines;  this  may 
be  unfavorable  to  the  performance  of  the  MMM. 

To  summarize  the  results,  Figures  1  through  4  present  the  min  and  max  Kendall-Tau 
test  p-value  for  each  replication  for  the  method  that  gives  the  best  fit. 

From  these  figures  and  specially  from  the  max  p-values,  it  may  be  deduced  that,  for  this 
data,  in  general  the  larger  the  max  p-value  (more  independence)  the  better  the  fit  with  the  MLE 
method.  Conversely  the  smaller  the  max  p-value  (less  independence)  the  better  the  fit  with  the 
MMM  method. 
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Table  7.  Correlation  -  20  Engines 


Correlation  Coefficients  for  20  enqines 

repl 

Failurel 
Failure2 

Failure3 

Failurel 

Failure? 

Failure3 

rep  16 

Failurel 
Failure? 
Failure^ 

Failurel 

Failure2 

Failure) 

1 

1 

-0  099201099 

1 

0  381725000 

1 

-0  325682845 

0 177625785 

1 

0  095575000 

-0  033370000 

1 

rep2 

Failurel 
Failure2 
Failure3 

Failurel 

Failure2 

Failure3 

rep  17 

Failurel 

Failure? 
Failure3 

Failurel 

Failure2 

Failure3 

1 

1 

-0  251973035 

1 

-0  214350000 

1 

0  460278255 

-0.055862733 

1 

-0  567440000 

0066732000 

1 

rep3 

Failurel 
Failure2 

Failure3 

Failurel 

Failure2 

Failures 

rep18 

Failurel 

Failure? 
Failure3 

Failurel 

Failure2 

Failures 

1 

1 

-0  292364468 

1 

-0  019240000 

1 

-0  171745200 

-0  161919234 

1 

0  089556000 

0  366561000 

1 

rep4 

Failurel 
Failure2 

Failure3 

Failurel 

Failure? 

Failure3 

rep  19 

Failurel 

Failure? 
Failure3 

Failurel 

Failure2 

Failures 

1 

1 

-0  266964793 

1 

0  062676000 

1 

0  375627368 

-0  020033442 

1 

-0  211740000 

-0  000530000 

1 

rep5 

Failurel 
Failure2 
Failure3 

Failurel 

Failure2 

Failure3 

rep20 

Failurel 

Failure? 
Farlure3 

Failurel 

Failure? 

Failure3 

1 

1 

-0  125786346 

1 

-0  227550000 

1 

-0  240386929 

-0  146477812 

1 

0225389000 

-0.371290000 

1 

rep6 

Failurel 

Failure? 
Failure3 

Failurel 

Failure2 

Failure3 

rep21 

Failurel 

Failure? 
Farlure3 

Failurel 

Failure? 

Failure3 

1 

1 

-0  049045517 

1 

-0  267500000 

1 

0 141463053 

0  466711164 

1 

0  140151000 

-0  088550000 

1 

rep7 

Failurel 

Failure  2 
Failure3 

Failurel 

Failure  2 

Failure3 

rep22 

Failurel 

Failure? 

Failure?. 

Failurel 

Failure? 

Failures 

1 

1 

-0  334293594 

1 

-0  237640000 

1 

-0  102335000 

0 108178021 

1 

0.051596000 

-0  117030000 

1 

rep8 

Failurel 

Failure  2 
Failure3 

Failurel 

Failure? 

Failure3 

rep23 

Failurel 

Failure? 
Failure3 

Failurel 

Failure2 

Failure3 

1 

1 

-0  362844249 

1 

-0  300890000 

1 

-0.153316474 

-0  014631447 

1 

0  077012000 

-0  043760000 

1 

rep9 

Failurel 

Failure? 
Failure3 

Failurel 

Failure2 

Failure3 

rep24 

Failurel 

Failure2 

Failures 

Failurel 

Failure? 

Failures 

1 

1 

0  138519930 

1 

-0  159170000 

1 

-0  383549222 

0  096601004 

1 

0  008203000 

0  087306000 

1 

rep  10 

Failurel 

Failure? 
Fatlure3 

Failurel 

Failure2 

Fajlure3 

rep25 

Failurel 
Failure2 
Failure3 

Failurel 

Failure2 

Failures 

1 

1 

-0  205353856 

1 

0  368302000 

1 

-0  202946234 

-0  062062633 

1 

-0  205760000 

-0  176000000 

rep  11 

Failurel 
Failure2 
Failures 

Failurel 

Failure2 

Failure3 

rep26 

Failurel 

Failure? 
Failures 

Failurel 

Failure2 

Fai!ure3 

1 

1 

-0  23980229 

1 

-0076760000 

1 

0  164227499 

-0  089182145 

1 

0025953000 

-0  17618 

1 

rep  12 

Failurel 

Failure? 
Fajlure3 

Failurel 

Failure2 

Failures 

rep27 

Failurel 
Failure2 

Failures 

Failurel 

Failure2 

Failure3 

1 

1 

-0  291025767 

1 

-0  059610000 

1 

-0283057458 

-0  216673298 

1 

0009859000 

0  421248000 

1 

rep13 

Failurel 

Failure? 
Failure3 

Failurel 

Failure2 

Failure3 

rep28 

Failurel 
Failure2 

Failure3 

Failurel 

Failure2 

Failure3 

1 

1 

0  010730000 

1 

-0  166470000 

1 

-0  148510000 

0  018339000 

1 

-0  408460000 

-0  182890000 

1 

rep  14 

Failurel 

Failure? 
Failure3 

Failurel 

Failure? 

Failure3 

rep29 

Failurel 

Failure? 
Failure3 

Failurel 

Failure? 

Failure3 

1 

1 

-0  431010000 

1 

-0  183030000 

1 

-0  284860000 

-0  144210000 

1 

-0  007070000 

-0  235020000 

1 

repIS 

Failurel 

Failure? 
Failure3 

Failurel 

Failure2 

Failure3 

rep30 

Failurel 
Failure2 

Failures 

Failurel 

Failure? 

Failure3 

1 

1 

-0  075920000 

1 

0011941000 

1 

0 111095000 

-0  381490000 

1 

0 144591000 

0 137677000 

1 
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Table  8.  Correlation  -  100  Engines 


Correlation  Coefficients  for  100  enqines 

repl 

Failir  e1 
Failire2 

F.nlif  e3 

Failir  e1 

Failire2 

Failir  e3 

rep16 

Failirel 
Failire2 
Failire3 

Failirel 

Failir  e2 

Failir  e3 

1 

1 

0  012345873 

1 

-0  343673022 

1 

-0  151985402 

0  027275574 

1 

-0  065777779 

-0  047403929 

1 

rep2 

Failirel 

Failure2 
Failure3 

Failirel 

Failire2 

Failir  e3 

rep  17 

Failirel 
Failure2 

Failir  e3 

Failirel 

Failire2 

Failir  e3 

1 

1 

■0  053286726 

1 

-0  095841896 

1 

-0  155893588 

-0  044274299 

1 

0 153445659 

-0  210881943 

r»p3 

Faili»*1 
Failtve? 
Failure3 

Failurel 

Failire2 

Failir  e3 

rep18 

Failirel 
Failir  e2 
Failir  e3 

Failirel 

Failure2 

Failir  e3 

1 

1 

-0.250749723 

1 

0052473163 

1 

0  062100274 

-0  094868699 

1 

0  041985512 

0  022318537 

1 

rep4 

Failirel 
Failure2 

Failure3 

Failirel 

Failire2 

Failir  e3 

rep19 

Failirel 
Failir  e2 
Failir  e3 

Failirel 

Failir  e2 

Failir  e3 

1 

1 

-0218109532 

1 

-0  062164561 

1 

-0  144770967 

-0  109527904 

1 

-0  133619298 

-0  160869201 

1 

repS 

Failirel 

Failure? 
Failir  e3 

Failirel 

Fallire2 

Failir  e3 

rep20 

Failirel 
Failir  e2 
Failir  e3 

Failirel 

Failure2 

Failir  e3 

1 

1 

-0  150008998 

1 

-0  083449627 

0 132158154 

-0  078813135 

1 

-0  003128899 

0  055670639 

1 

rep6 

Failirel 
Failir  e2 
Failure3 

Failirel 

Failtre2 

Faituro3 

rep21 

Failirel 
Failir  e2 
Failir  e3 

Failirel 

Failir  e2 

Failir  e3 

1 

1 

-0079147073 

1 

-0  110695458 

1 

-0  240225116 

-0  170526078 

1 

-0  170345154 

-0  103486186 

1 

rep7 

Failirel 
Failire2 

Failure3 

Failirel 

Fallire2 

Failir  e3 

rep22 

Failirel 
Failure2 
Failir  e3 

Failirel 

Failir  e2 

Failir  e3 

1 

1 

-0  065957419 

1 

-0  180701469 

1 

-0  065309916 

-0  149967098 

1 

0 156143924 

-0  219732137 

l 

rap8 

Failirel 
Fajtire2 

Failure3 

Failurel 

Failir  e2 

Failir  e3 

rep23 

Failirel 

Failure? 
Failir  e3 

Failirel 

Failir  e2 

Failir  e3 

1 

1 

-0  087122634 

1 

-0073156230 

1 

0  009267376 

-0  198193315 

1 

-0  171758743 

-0  083586000 

1 

rep9 

Failirel 
Failire2 
Failire3 

Failirel 

Failir  e2 

Failir  e3 

rep24 

Failirel 
Failir  e2 
Failir  e3 

Failirel 

Failir  e2 

Failir  e3 

1 

1 

-0  024708152 

1 

-0  063590173 

1 

-0  088238686 

-0  153996216 

1 

-0  106654719 

-0  162453112 

1 

rep  10 

Failirel 
Failire2 
Failure3 

Failirel 

Failire2 

Failir  e3 

rep25 

Failirel 
Failir  e2 

Failir  e3 

Failirel 

Failir  e2 

Failir  e3 

1 

1 

-0  085599667 

1 

-0  084052612 

1 

-0  208242972 

-0  04475226 

1 

0049055891 

-0  063152872 

1 

rap11 

Failirel 
Failire2 
Failure3 

Failirel 

Failir  e2 

Failir  e3 

rep26 

Failirel 
Failir  e2 
Failir  e3 

Failirel 

Failire2 

Faitire3 

1 

1 

-0  030401682 

1 

-0  066404142 

1 

0  019141802 

0  012598681 

1 

0  099904638 

-0219739207 

1 

rep12 

Failirel 
Failire2 

Failure3 

Failirel 

Failure2 

Failures 

rep27 

Failirel 
Failir  e2 
Failir  e3 

Failirel 

Failir  e2 

Failir  e3 

1 

1 

-0  243204236 

1 

-0  056026419 

1 

0  0168114:7 

-0  06918138 

1 

-0  077399648 

-0  002855117 

1 

rep  13 

Failirel 
Failire2 
Failure3 

Failirel 

Failir  e2 

Failir  e3 

rep28 

Failirel 
Failire2 

Failir  e3 

Failirel 

Failir  e2 

Failir  e3 

1 

1 

-0250749723 

1 

-0  081782938 

1 

0  062100274 

-0  094868699 

1 

-0  061450025 

0  020642837 

1 

rep  14 

Failirel 
Failure2 

Failire3 

Failirel 

Failire2 

Failir  e3 

rop29 

Failirel 

Failure? 
Failir  e3 

Failirel 

Failir  e2 

Failir  e3 

1 

1 

-0  064039892 

1 

-0  243170112 

1 

-0  077445912 

-0  077941408 

1 

0  093441256 

-0  10701961 

1 

repIS 

Failirel 
Failire2 
Faflure3 

Failirel 

Failir  e2 

Failir  e3 

rep30 

Failirel 

Failire? 
Failire3 

Failirel 

Failire2 

Failir  e3 

1 

1 

-0  098178961 

1 

-0  23939123 

1 

0  077241660 

-0  170746677 

1 

-0  020603539 

-0  257213197 

1 

22 


Table  9.  Kendall  -Tau  Test  -  20  Engines 
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KENDALL-TAU  TEST  OF  INDEPENDENCE  FOR  THE  SET  OF  20  ENGINES 

tau  12 

pvalue 

tau13 

pvalue 

tau23 

pvalue 

mean 

min 

max 

best  method 

1 

-0  084 

0.604 

-0.242 

0.136 

0.126 

0.436 

0.392 

0.136 

0.604 

m.m.m 

2 

-0.211 

0.194 

0.221 

0.173 

-0.062 

0.697 

0.355 

0.173 

0.697 

m  m.m 

3 

-0  116 

0.475 

-0.158 

0.330 

-0.095 

0.559 

0.455 

0.330 

0.559 

m.m.m 

4 

-0.179 

0.270 

0.105 

0.516 

0.126 

0.436 

0.408 

0.270 

0.516 

m.l.e. 

5 

-0.126 

0.436 

-0.295 

0.069 

-0.074 

0.650 

0.385 

0.069 

0.650 

m.m.m 

6 

0.105 

0.516 

0.105 

0.516 

0.221 

0.173 

0.402 

0.173 

0.516 

m.m.m 

7 

-0.232 

0.153 

-0.011 

0.948 

0.084 

0.604 

0.568 

0.153 

0.948 

m.l.e. 

8 

-0.221 

0.173 

-0.063 

0.697 

-0.168 

0.299 

0.390 

0.173 

0.697 

m.l.e. 

9 

0.063 

0.697 

-0.295 

0.069 

0.179 

0.270 

0.345 

0.069 

0.697 

m.l.e. 

10 

-0.011 

0.948 

-0.274 

0.092 

-0.126 

0.436 

0.492 

0.092 

0.948 

m.m.m 

11 

-0.147 

0.364 

0.010 

0.948 

-0.063 

0.697 

0.670 

0.364 

0.948 

m.l.e. 

12 

-0.116 

0.475 

-0.158 

0.330 

-0.095 

0.559 

0.455 

0.330 

0.559 

m.l.e. 

13 

-0.042 

0.795 

-0.179 

0.270 

0.147 

0.364 

0.476 

0.270 

0.795 

m.l.e. 

14 

-0.263 

0.105 

-0.179 

0.270 

-0.179 

0.270 

0.215 

0.105 

0.270 

m.m.m 

15 

-0.168 

0.299 

0.232 

0.153 

-0.326 

0.044 

0.166 

0.044 

0.299 

m.m.m 

16 

0.189 

0.243 

0.221 

0.173 

-0.021 

0.897 

0.438 

0.173 

0.897 

m.l.e. 

17 

-0.126 

0.436 

-0.400 

0.014 

0.116 

0.475 

0.308 

0.014 

0.475 

m.l.e 

18 

0.105 

0.516 

0.053 

0.746 

0.105 

0.746 

0.669 

0.516 

0.746 

m.l.e. 

19 

-0.126 

0.436 

-0.063 

0.697 

0.011 

0.948 

0.694 

0.436 

0.948 

m.m.m 

20 

-0.116 

0.475 

0.221 

0.173 

-0.200 

0.218 

0.289 

0.173 

0.475 

m.l.e. 

21 

-0.168 

0.218 

0.274 

0.092 

-0.095 

0.559 

0.289 

0.092 

0.559 

m.m.m 

22 

-0.211 

0.194 

-0.095 

0.559 

0.042 

0.795 

0.516 

0.194 

0.795 

m.l.e. 

23 

-0.189 

0.243 

0.242 

0.136 

-0.042 

0.795 

0.391 

0.136 

0.795 

m.l.e. 

24 

-0.074 

0.650 

-0.116 

0.475 

0.032 

0.846 

0.657 

0.475 

0.846 

m.l.e. 

25 

0.147 

0.364 

-0.253 

0.119 

-0.053 

0.746 

0.410 

0.119 

0.746 

m.l.e. 

26 

-0.042 

0.795 

0.021 

0.897 

-0.053 

0.746 

0.813 

0.746 

0.897 

m.m.m 

27 

-0.074 

0.650 

0.147 

0.364 

0.105 

0.516 

0.510 

0.364 

0.650 

m.m.m 

23 

-0.042 

0.795 

-0.263 

0.105 

-0.084 

0.604 

0.501 

0.105 

0.795 

m.l.e. 

29 

-0.074 

0.650 

-0.105 

0.516 

-0.105 

0.516 

0.561 

0.516 

0.650 

m.l.e. 

30 

-0  095 

0.559 

0.032 

0.846 

0032 

0.846 

0.750 

0.559 

0  846 

m.l.e. 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

-0.08804 

045767 

-0.04211 

0.38100 

-0.01714 

0.55823 

0.46563 

0.24566 

0.69414 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

0.11499 

0.22407 

0.19334 

0.29522 

0 12607 

023159 

0  15545 

0.18088 

0  18533 

23 


Table  10.  Kendall  -Tau  Test  - 100  Engines 


KENDALL-TAU  TEST  OF  INDEPENDENCE  FOR  THE  SET  OF  100  ENGINES 

tau  12 

pvalue 

tau13 

pvalue 

tau23 

pvalue 

mean 

min 

max 

best  method 

1 

0.029 

0.668 

-0.136 

0.045 

-0.072 

0.292 

0.335 

0.045 

0.668 

m.l.e 

2 

-0.093 

0.173 

-0.181 

0.008 

0.003 

0.967 

0.382 

0.008 

0.967 

m.l.e. 

3 

-0.181 

0.008 

0.003 

0.967 

-0.025 

0.708 

0.561 

0.008 

0.967 

m.l.e. 

4 

-0.144 

0.034 

-0.792 

0.428 

-0.064 

0.347 

0.270 

0.034 

0.428 

m  m.m 

5 

-0.089 

0.190 

0.088 

0.192 

-0.036 

0.592 

0.325 

0.190 

0.592 

m.m.m 

6 

-0.138 

0.042 

-0.166 

0.015 

-0.098 

0.150 

0.069 

0.015 

0.150 

m.m.m 

7 

-0.028 

0.681 

-0.059 

0.381 

-0.064 

0.344 

0.469 

0.344 

0.681 

m.l.e. 

8 

-0.036 

0.592 

-0.053 

0.435 

-0.145 

0.033 

0.353 

0.033 

0.592 

m.m.m 

9 

-0.025 

0.708 

-0  009 

0.896 

-0.109 

0.108 

0.570 

0.108 

0.896 

m.l.e. 

10 

-0.026 

0.703 

-0.110 

0.104 

-0.019 

0.775 

0.527 

0.104 

0.775 

m.m.m 

£ 

11 

-0.002 

0.976 

0.042 

0.536 

0.026 

0.703 

0.738 

0.536 

0.976 

m  I.e. 

3 

12 

-0.132 

0.052 

0.064 

0.347 

-0.017 

0.798 

0.399 

0.052 

0.798 

m.l.e. 

a 

13 

-0.181 

0.008 

0.003 

0.967 

-0.025 

0.708 

0.561 

0.008 

0.967 

m  I.e. 

c 

88 

14 

-0.081 

0.231 

-0.069 

0.308 

-0.018 

0.789 

0443 

0.231 

0.789 

m.l.e. 

_u 

15 

-0.051 

0.449 

0.095 

0.162 

-0.093 

0.173 

0.261 

0.162 

0.449 

m.m.m 

a 

16 

-0.250 

0.000 

-0.060 

0.378 

-0.039 

0.568 

0.315 

0.000 

0.568 

m  I.e. 

a: 

17 

-0.070 

0.303 

0.011 

0.868 

-0.158 

0.020 

0.397 

0.020 

0.868 

m.l.e. 

IB 

-0.059 

0.385 

0.016 

0.816 

-0.004 

0.953 

0.718 

0.385 

0.953 

m.m.m 

19 

-0.074 

0.273 

-0.053 

0.435 

-0.134 

0.048 

0.252 

0.048 

0.435 

m.l.e. 

20 

-0.055 

0.415 

-0.022 

0.748 

-0.023 

0.739 

0.634 

0.415 

0.748 

m  I.e. 

21 

-0.058 

0.394 

-0.120 

0.077 

-0.068 

0.314 

0.262 

0.077 

0.394 

m.m.m 

22 

-0.104 

0.126 

0.122 

0.072 

-0.198 

0  004 

0.067 

0.004 

0.126 

m.l.e. 

23 

-0.046 

0.501 

-0.161 

0.018 

-0.135 

0.046 

0.188 

0.018 

0.501 

m.l.e. 

24 

-0.093 

0.173 

-0.118 

0.082 

-0.058 

0.391 

0.215 

0.082 

0.391 

m.l.e. 

25 

0.022 

0.748 

0.071 

0.297 

-0.039 

0.568 

0.538 

0.297 

0.748 

m.l.e. 

26 

-0.031 

0.651 

0.041 

0.544 

-0.210 

0.002 

0.399 

0.002 

0.651 

m.l.e. 

27 

0.006 

0.929 

-0.051 

0.457 

-0.002 

0.976 

0.787 

0.457 

0.976 

m.l.e. 

28 

-0.029 

0.668 

-0.005 

0.938 

0.017 

0.803 

0.803 

0.668 

0.938 

m.m.m 

29 

-0.160 

0.019 

-0.022 

0.743 

-0.020 

0.770 

0.511 

0.019 

0.770 

m  I.e. 

30 

-0.189 

0.005 

0.002 

0.981 

-0.154 

0.024 

0.337 

0.005 

0.981 

mm  m 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

-007884 

037010 

-0.05429 

0.44144 

-0  06601 

045688 

0.42281 

0.14569 

069136 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

0  06784 

0  30264 

0 16079 

033413 

006379 

033785 

0 19490 

0 18499 

0.24739 

24 
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Figure  1.  Maximum  p- Value  for  20  Engines 
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Figure  2.  Minimum  p- Value  for  20  Engines 
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Figure  3.  Maximum  p- Value  for  100  Engines 
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Figure  4.  Minimum  p- Value  for  100  Engines 
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V.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


A.         SUMMARY 

The  analytical  approach  taken  to  develop  this  thesis  and  topics  covered  has  been  the 
following: 

1.  Develop  a  program  in  S-Plus  to  generate  the  failure  times  for  an  engine  with  3 
subsystems,  where  the  distribution  of  the  first  failure  times  of  these  subsystems  is 
Weibull  and  the  engine  is  not  restored  to  new  condition  after  each  engine  failure. 

2.  Estimate  the  parameters  for  the  Weibull  distribution  that  fits  the  generated  engine 
failure  times  with  the  Maximum  Likelihood  Estimator  and  a  Modified  Method  of 
the  Moments  using  commercial  spread-sheet  software. 

3 .  Measure  the  accuracy  of  these  two  methods  with  the  Kolmogorov-Smirnov  Test. 

4.  Apply  the  same  analysis  to  a  given  set  of  real  engine  failure  times. 

5.  Compare  the  accuracy  of  the  two  methods,  and  discuss  reasons  for  their 
differences. 


B.         CONCLUSIONS 

The  two  methods  in  general  give  good  results  when  they  are  applied  to  fit  a  Weibull 
distribution  to  the  failure  times  generated  in  the  simulation.  The  resulting  fit  for  these  failure 
times  was  better  with  the  Maximum  Likelihood  Estimator  more  often  than  with  the  Modified 
Moment  Method.  For  this  reason,  and  the  fact  that  the  MLE  method  has  superior  theoretical 
properties  in  the  presence  of  independence,  the  author  recommends  the  MLE  method  to 
estimate  the  parameters  for  scenarios  similar  to  those  simulated  in  this  Thesis. 

With  these  conclusions  it  is  important  to  point  out  the  following  considerations: 
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The  MMM  has  the  advantage  that  fewer  calculations  are  necessary  to  get  the  results. 
In  the  case  of  the  Macros  developed  in  this  study,  only  the  initial  values  and  two  columns  of 
calculations  are  necessary.  Four  columns  of  calculations  are  necessary  in  the  MLE  case. 

There  was  a  wide  range  of  values  for  the  parameters  of  the  fitted  distributions  across  all 
the  replications.  This  is  displayed  in  Table  1 1 . 

Table  11.  Parameters  Range 


20  Engines 

MLE 

MMM 

0.961 1<=  shape  <=  1.5238 

0.9857  <=  shape  <=  1.4616 

223.28  <=  scale  <=  330.53 

236.97  <=  scale  <=  329.43 

100  Engines 

MLE 

MMM 

1.1002  <=  shape  <=  1,3422 

1.1737  <=  shape  <=  1.4485 

255.64  <=  scale  <=  301.25 

254.07  <=  scale  <=  305.76 

The  dependence  obtained  in  the  simulation  is  not  as  significant  as  was  expected. 

C.         RECOMMENDATIONS 

These  two  methods  should  be  studied  when  greater  dependence  exists  among  the  data 
and  when  the  dependency  is  introduced  in  a  different  manner  than  used  in  this  thesis. 
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APPENDIX  A.  PROGRAM  TO  GENERATE  FAILURE  TIMES 


fuiiction(n) 

{ 

#  generate 

#  This  function  generates  times  randomly  in  a  manner  that  resembles  the  way  in  which 

#  engines  failures  occur  in  the  real  world. 

totalreplic  <-  NULL 
for(count  in  l:n)  { 

#  Step  1 

#  Three  Weibull  random  numbers  are  generated,  each  one  represents  the  first  failure  time 

#  that  each  one  of  the  components  could  have.  These  values  constitute  the  vector  faultsl. 

#  Three  Weibull  random  numbers  are  generated,  each  one  represents  the  second  failure  time 

#  that  each  one  of  the  components  could  have.  These  values  constitute  the  vector  faults2. 

#  Three  Weibull  random  numbers  are  generated,  each  one  represents  the  third  failure  time 

#  that  each  one  of  the  components  could  have.  These  values  constitute  the  vector  faults3. 

xll<-rweibull(l,  1.2,900) 

xl2<-rweibull(l,  1.3,800) 

xl3<-rweibull(l,  1.4,600) 

faultsl  <-c(xll,  xl2,  xl3) 

x21<-rweibull(l,  1.2,900) 

x22<-rweibull(l,  1.3,800) 

x23  <-  rweibull(l,  1.4,600) 

faults2  <-  c(x21,  x22,  x23) 

x3K-rweibull(l,  1.2,900) 

x32  <- rweibull(l,  1.3,800) 

x33<-rweibull(l,  1.4,600) 

faults3<-c(x31,  x32,  x33) 

failures  <-  matrix(c(faultsl.  faults2,  faults3),  ncol  =  3) 

#  Step  2 

#  The  minimum  value  of  the  vector  faultsl  is  selected,  this  value  is  the  time  of  the  first 

#  failure  of  the  engine.  The  program  creates  a  new  vector  newfail2  that  is  the  difference 

#  between  the  vector  faultsl  and  the  minimum  value  selected,  it  is  clear  that  this  new  vector 

#  has  a  component  that  is  0,  this  component  is  replaced  in  the  vector  newfail2  for  the 

#  element  of  the  vector  faults2  that  match  with  the  component  which  had  the  first  failure. 

#  Choose  the  first  failure  time,  that  is  the  minimum 

#  time  from  the  vector  faultsl 

fl  <-  min(faultsl) 

#  the  next  step  is  to  choose  the  second  failure  time,  to  make  this  first 

#  the  above  value  min(faultsl)  must  be  substituted  for  a  new  time 

#  generated  in  the  vector  faults2  and  precisely  in  the  same  component. 

#  After  this  the  new  minimum  is  selected 

newfail2<- (faultsl  -fl) 
ii  <-  match(0.  newfail2) 
newfail2[ii]  <-  faults2[ii] 

#  Step  3 

#  At  this  point  the  engine  has  two  components  that  have  accumulated  working  hours  and 

#  one  that  has  been  just  repaired.  The  vector  newfail2  has  two  elements  that  have  the 
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#  remainder  of  the  theoretical  first  failure  for  these  components  and  one  element  that  is  the 

#  theoretical  second  failure  for  the  other  component.  So,  again,  the  minimum  element  of 

#  this  vector  is  selected  and  it  is  the  second  failure  time  of  the  engine. 

f2  <-  min(newfail2) 

#  Step  4 

#  The  program  creates  a  new  vector  called  newfaiO  that  is  the  difference  between  the  vector 

#  newfail2  and  this  minimum  value  selected,  it  is  clear  that  this  vector  has  a  component  that 

#isO. 

#  Choose  the  third  failure  time,  here  there  are  two  possible 

#  cases,  that  the  second  failure  is  in  the  same  component  that 

#  the  first  one  or  not,  these  two  condition  are  programmed  in  this  order 
# 

newfail3  <-  (newfail2  -  f2) 
iii  <-  match(0,  newfail3) 
if(ii  =  iii)  { 

newfail3[iii)  <-faults3[iii] 

D  <-  min(newfail3) 


: 


else  if(ii  !=  iii)  { 
newfail3[iii]  <-faults2[iii] 

f3  <-  min(newfail3) 

enginefails  <-  c(fl,  f2,  D) 
#cumenginefails  <-  cumsum(enginefails) 

#engfails  <-  matrix(c(failures,  enginefails,  cumenginefails),  byrow  = 
#   T,ncol  =  3) 

totalreplic  <-  rbind(totalreplic,  (enginefails)) 

totalreplic 
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Row  Names 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 


V1 

V2 

V3 

88.980852932 

7.374101636 

166.612433434 

57.770945753 

452.164077047 

671.533220087 

348.363054836 

13.725217444 

33.719418711 

382.257870284 

168.795667771 

119.352260671 

57.145359689 

117.795729216 

228.544288266 

251.488400993 

417.099715694 

206.169150245 

703.336492627 

16.866142294 

168.320579106 

85.174113550 

960.378668893 

561 .250692879 

126.379607253 

166.962891751 

410.881188239 

9.871149363 

30.987527462 

1171.995737937 

184.016586935 

55.183363300 

170.715797827 

19.085994032 

206.177991742 

124.656084035 

171.708369596 

7.411273439 

458.366290751 

957.877021060 

334.941605556 

252.553811069 

88.483061291 

129.532846338 

55.052226186 

97.662694658 

64.764644477 

521.223680534 

213.546020387 

8.955316669 

314.364020817 

420.340591055 

182.036946622 

345.506361499 

46.823610226 

411.857588415 

109.450083602 

25.300667331 

361.149199434 

744.757337329 

Column  1            Column  2             Column  3 

Column  1 
Column  2 
Column  3 

1 

-0.099201099 

1 

-0.325682845 

0.177625785 

1 
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Row  Names 
1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 


V1 

V2 

V3 

450.634670353 

134.620284266 

763.602013650 

395.264314931 

300.696009190 

201 .433889226 

467.906036969 

301.510351673 

22.166715862 

200.004796884 

261.419709923 

234.019837022 

62.187250449 

351.566316602 

293.954533661 

188.420555520 

530.441439163 

446.499347177 

517.045829232 

212.650323103 

4.37229061 2 

104.406621033 

91  440692508 

286.231002289 

185  580658025 

189.658931128 

461 .45231 7203 

278.805459607 

142.687073420 

63.785901768 

242.374076871 

140.003133652 

57.571537504 

191.638384699 

35.982496826 

135.685578893 

115.651429115 

807.515319374 

56.217725189 

318.738490247 

220.511514630 

96.958868093 

389.324028552 

121.699393963 

391 .967785281 

108.549643630 

833.642731237 

587.010089214 

647.283140422 

270.006356536 

600.593877204 

135.277879814 

235.260576261 

294.610909130 

167.565552866 

1 9.382052477 

80.687152199 

620.417102151 

43.940840941 

348.339229735 

549.652520297 

147.263663645 

402.329146256 

160.858623029 

198.697751449 

131  936073250 

269.748636292 

216.407382425 

272.564254146 

119.882883842 

42.753769287 

568  105938757 

290.367131217 

46.018913125 

470.135356798 

111.526180647 

59.040162467 

442.572648090 

315.869668084 

112.068427109 

799073849390 

258.473619251 

740.012832046 

730.391096774 

936.074568617 

525.344669063 

53.458836253 

38.275293391 

264  758222487 

109.033459770 

200.399299969 

130.448227050 

564.453719467 

1 6.387049764 

689.920420726 

124.595174314 

282.346963237 

302.451537677 

17.105846707 

298.005804703 

792.704752916 

277.592370478 

488.070794225 

600.835486551 

77.330402620 

651.631129729 

186.384979015 

302.942800237 

488.356688854 

282.631614706 

1 2.094239822 

251.322489919 

639.047414057 

315  553538531 

102.245119335 

121.065856465 

248  600567386 

120.292743494 

281.418670027 

6.51  2039507 

394.628880256 

329.862437764 

42 172683323 

275.823361976 

400.377074953 

205.279868273 

89.822151642 

568.897212277 

126.945901458 

413.319725633 

280.829982419 

74.621453179 

172.347938766 

138.081515648 

169  690890193 

295.648493372 

49.777202258 

271  146965769 

427  883569557 

120.358845277 

86  140070778 

106.453003164 

336.696295308 

392  945334315 

439.747178352 

145.848786206 

357  670884257 

95424425759 

158.394437081 

233  214847248 
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Row  Names 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 


V1 

V2 

V3 

78  318115166 

233  292392917 

150.556841342 

400.322169035 

703.176710554 

131  985760999 

284.082094482 

700.398323845 

54.502213913 

673.813459742 

628.345709323 

27.810511050 

416.157442911 

309.802781935 

61.701311275 

502.727434228 

208.984463635 

7  000471997 

395.249167889 

181.615779362 

246.169600960 

532.055233005 

282.774380172 

51  530633142 

256.526079351 

265.073575146 

366.480188363 

791.013013348 

483.679778572 

383.122947523 

116.674678096 

203.429431816 

501 .860086251 

216.019262031 

50.783576381 

346.268316228 

652.055312871 

297.710549841 

19.143684755 

57.143580125 

444.863422850 

53.260846841 

305.430041853 

64.522458661 

344.328881540 

134.875969833 

153  848996778 

74.573546508 

207.795329813 

197.474623424 

36.572283433 

55.079317773 

148.585617029 

153.485311966 

332.356937878 

183.167123218 

2.771444009 

33.292174799 

281.412899135 

521  055842681 

329.237256924 

624.659923834 

410.145804201 

134.999172015 

924.566774196 

959.801158180 

310.283780658 

196.949468983 

61 .420855422 

145.927518733 

98.091275956 

105.358948583 

96.240971550 

448.577731504 

68.34921 2862 

504.457185330 

3.995392340 

54.819812986 

779.185987694 

63.103630524 

510.654006981 

64.295309451 

227.558484641 

681.717084473 

132.919226467 

80.090112060 

305.523448635 

225.006666940 

354.405022103 

115.791651829 

322.833702315 

15.741543588 

273.716914308 

552.256221658 

238.880493230 

103.236641646 

234.905588128 

267.774447132 

384.884639091 

241.821803683 

27.660947307 

35.559153707 

360.121368725 

54.4003091 93 

74.75701 2406 

109.375105806 

169.046509653 

306.549361255 

211.877188712 

54.427577933 

47.413167726 

328.701690509 

233.519929404 

411.204384691 

271.788277670 

132.777355320 

359.728088427 

54.715071977 

52.466467632 

253.216046834 

157.314280817 

49.880943701 

746.100150202 

93.608619218 

110.699690363 

573390687287 

218.283430910 

434.282741182 

109.443120744 

1106.554817220 

178.424617471 

107.752314459 

179  703059410 

96.191152623 

63  123472845 

220.753490182 

358.476353401 

438.979104011 

97.987549536 

57.556717482 

152.972644520 

126.992796611 

173.350825362 

1004.747325381 

656.574251587 

250.152057881 

312.989740534 

189  330584986 

142.057479837 

63.629475350 

Column  1 

Column  2 

Column  3 

Column  1 
Column  2 
Column  3 

1 

0.012345873 

1 

-0  151985402 

0  027275574 

1 
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APPENDIX  B.  MACRO  20  ENGINES 


MACRO  FOR  20  ENGINES 

thesis  Macro 

Macro  recorded  bv  FERNANDO  GU1MIL 


Sub  thesisO 
Sheets  Add 

Sheets("Sheet2").Select 
Range("Al:C2").Select 

ActiveCell.FormulaRlCl  =  "MAXIMUM  LIKELIHOOD  ESTIMATOR" 
Range("Al:C2").Select 
Range("A2")Activate 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAJignment  =  xlCenter 

Vertical  Ahgnment  =  xlBottom 

.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
With  Selection 

.HorizontalAlignment  =  xlCenterAcrossSelection 

.Vertical  Alignment  =  xlBottom 

WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 
Range("Dl").  Select 

ActiveWindow  SmallScroll  ToRight:=10 
With  Toolbars(8) 

.Left  =  241 

Top  =  228 
End  With 

ActiveWindow.LargeScrollToRight:=-l 
Range("Pl:Rl").Select 

ActiveCell.FormulaRlCl  =  "MODIFIED  MOMENT  METHOD" 
Range("Pl:Rl").Select 
Range("Ql").Activate 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

.Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
With  Selection 

HorizontalAlignment  =  xlCenterAcrossSelection 

VerticalAlignment  =  xlBottom 
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.WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 

Range("Nl").Select 

ActiveWindow.LargeScroll  ToRight:=-l 
With  TooIbars(8) 

.Left  =  436 

.Top  =  226 
End  With 

Range("A2").Select 

ActiveCell.FormulaRlCl  =  "shape  parameter" 
With  ActiveCell.Characters(Start:=l,  Length: =15).Font 

.Name  =  "Arial" 

FontStyle  =  "Bold" 

.Size  =  10 

Strikethrough  =  False 

Superscript  =  False 

Subscript  =  False 

OutlineFont  =  False 

Shadow  =  False 

Underline  =  xlNone 

Colorlndex  =  xl  Automatic 
End  With 
Range("B3").  Select 

ActiveCell.FormulaRlCl  =  "initial  beta" 
Range("B3").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
Range("B5").  Select 

ActiveCell.FormulaRlCl  =  "final  beta" 
Range("B5").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("A7:C7").Select 
ActiveCell.FormulaRlCl  =  "scale  parameter" 
Range("A7:C7").Select 
Rangc("B7").  Activate 
Selection  Font. Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 
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.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
Range("A7").Select 
ActiveCell.ForniulaRlCl  =  "" 
Range("B7").  Select 

ActiveCell.ForniulaRlCl  =  "scale  parameter" 
Range("B8").  Select 
ActiveCell.FormulaRlCl  =  "landa" 
Range("B8").  Select 
Selection.Font.Bold  =  True 
With  Selection 

Horizontal  Alignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Active  Window.  SmallScroll  ToRight:=0 
Range("A2:C2").Select 
Selection.  Copy 

Active  Window.  SmallScroll  ToRight:=l  1 
Range("Q2").Select 
ActiveSheet.Paste 
Application.  CutCopy  Mode  =  False 
Selection.  Clear 

ActiveWindow.LargeScroll  ToRight:=-2 
Range("A2:C2").Select 
Selection.  Copy 

ActiveWindow.  Small  Scroll  ToRight:=9 
Range("P2").  Select 
ActiveSheet.Paste 
Range("04").  Select 

ActiveWindow.LargeScroll  ToRight:=- 1 
Range("B3").  Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.  Small  Scroll  ToRight:=l  1 
Range("Q3").  Select 
ActiveSheet.Paste 

ActiveWindow.LargeScroll  ToRight  --2 
Range("B5").Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.  SmallScroll  ToRight:=12 
Range("Q5").  Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 
ActiveSheet.Paste 
Range("R7").  Select 

ActiveWindow  LargeScroll  ToRight:=-2 
Range("B5").Select 
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Application.  CutCopyMode  =  False 

Selection. Copy 

ActiveWindow.SmallScroll  ToRight:=12 

Range("Q5").  Select 

ActiveSheet.Paste 

ActiveWindow.LargeScroll  ToRight:=-2 

Range("A7:C8").Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=l  1 

Range("P7").  Select 

ActiveSheet.Paste 

Range("N7").  Select 

ActiveWindow.LargeScroll  ToRight:=-2 

Range("D  11").  Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

Range("D7").Select 

ActiveSheet.Paste 

Range("B8").  Select 

ActiveWindow.SmallScroll  ToRight:=13 

ActiveWindow.LargeScroll  ToRight:=-2 

ActiveWindow.SmallScroll  ToRight:=T4 

ActiveWindow.LargeScroll  ToRight~-2 

Range("El").Select 

Application.  CutCopyMode  =  False 

ActiveCell.FormulaRlCl  =  "Number  of  observations" 

Range("El").Select 

Selection.Font.Bold  =  True 

With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
With  Selection 

HorizontalAlignment  =  xlCenterAcrossSelection 

Vertical  Alignment  =  xlBottom 

WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 

Range("E2").Select 
AcuveCell.FormulaRlCl  =  "n" 
Range("E2").Select 
With  Toolbars(8) 

Left  =  536 

Top  -  23 
End  With 

Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 
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Vertical  Alignment  =  xlBottom 
WrapText  =  False 
Orientation  =  xlHorizontal 

End  With 

Sheets("Sheetl").  Select 

Range("B2:B21").  Select 

Selection.  Copy 

Sheets("Sheet2").  Select 

Range("D10").Select 

ActiveSheet.Paste 

Sheets("Sheetl").Select 

Range("C2:C21")-Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

Sheets("Sheet2").  Select 

RangefE  10").  Select 

ActiveSheet.Paste 

Sheets("Sheetl").  Select 

Range("D2:D21").  Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

Sheets("  Sheet2 ").  Select 

Range("F10").Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

Sheets("Sheetl").Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

Sheets("Sheet2").Select 

ActiveSheet.Paste 

Range("D10:D29").Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

Range("G10").Select 

ActiveSheet.Paste 

Range("E10:E29").Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

Active  Window.  SmallScroll  Down:=16 

Range("G30").Select 

ActiveSheet.Paste 

ActiveWindow.LargeScroll  Down:=-l 

Range("F10:F29").Select 

Application.CutCopyMode  =  False 

Selection.  Copy 

Active  Window.  SmallScroll  Down:=35 

Range("G50").  Select 

ActiveSheet.Paste 

Range("D42").  Select 

ActiveWindow.LargeScroll  Down:=-2 

Range("G10:G69").Select 
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Application.  CutCopyMode  =  False 

Selection.  Copy 

Range("H10").Select 

ActiveSheet.Paste 

Application.  CutCopyMode  =  False 

Selection. Sort  Keyl:=Range("H10"),  Order l:=xl Ascending,  Header: 

xlGuess,  OrderCustom:=l,  MatchCase:=False,  Orientation—  _ 

xlTopToBottom 
Range("B  10").  Select 
ActiveCell.FormulaRlCl  =  "theta" 
Range("B10").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
Range("E3").  Select 
ActiveCell.FormulaRlCl  =  "60" 
Range("E3").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  -  xlCenter 

VerticalAlignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Active  Window.  ScrollRow=  1 
RangefGl").  Select 

ActiveCell.FormulaRlCl  =  "Failures  times  Ti" 
RangefGl").  Select 
With  Selection 

HorizontalAlignment  =  xl  Center  AcrossSelection 

VerticalAlignment  =  xlBottom 

WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 

Selection.Font.Bold  =  True 
Range("G9").  Select 
ActiveCell.FormulaRlCl  =  "Ti" 
Range("G9").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

VerticalAlignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
Selection.  Copy 
Range("H9").  Select 
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ActiveSheet.  Paste 

Range("G9").Select 

Application.  CutCopyMode  =  False 

Selection.  Clear 

ActiveWindow.  Small  Scroll  ToRight:=3 

Range("I9").  Select 

ActiveCell.FormulaRlCl  =  "Ln  Ti" 

Range("I9").  Select 

Selection.Font.Bold  =  True 

With  Selection 

Horizontal  Alignment  =  xl  Center 

.Vertical  Alignment  =  xlBottom 

.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("I10").Select 

ActiveCell.FormulaRlCl  =  "=LN(RC[-l])n 
Range("I10").Select 
Selection. AutoFill  Destination:=Range("I10:I69"),  Type: 

xlFillDefault 
Range("I10:I69").Select 
Range("H70").  Select 
ActiveCell.FormulaRlCl  =  "SUM" 
Range("H70").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

.Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
Selection.  Copy 
Range("I70").  Select 
ActiveSheet.Paste 
Range("H71").  Select 
Application.  CutCopyMode  =  False 
ActiveCell.FormulaRlCl  =  "=SUM(R[-61]C:R[-2]Q" 
Range("I71").Select 

AcuveCell.FormulaRlCl  =  "=SUM(R[-61]C:R(-2]C)" 
Range("I72").Select 
ActiveWindow.LargeScroll  Down:=-3 
Range("J9").Select 
AcuveCell.FormulaRlCl  =  "TiAB4" 
Range("J9").Select 
Selecuon.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

VerticalAlignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 


43 


Active  Window.LargeScroll  ToRight:=- 1 

Active  Window.  Small  Scroll  ToRight:=2 

Range("J10").Select 

ActiveCell.FormulaRlCl  =  "=RC[-2]AR[-6]C[-8]" 

Range("J10").Select 

Selection. AutoFill  Destination:=Range("J10:J69"),  1ype:=  _ 

xlFillDefault 
Range("J10:J69").Select 
Selection.Replace  What:- 'AB*M,  Replacement:- ,AB4",  LookAt:=xlPart 

SearchOrder:=xlByRows,  MatchCase:=False 
Range("K9").  Select 

Active  Window.  SmallScroll  ToRight:=2 
ActiveCell.FormulaRlCl  =  "TiAB6" 
Range("K9").  Select 
Selecuon.Font.Bold  =  True 
With  Selection 

.HorizontalAhgnment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("K10").Select 

ActiveCell.FormulaRlCl  =  "=RC[-3]AR[-4]C[-9]M 
Range("K10").Select 
Selection. AutoFill  Destination:=Range("K10:K69"),  Type:=  _ 

xlFillDefault 
Range("K10:K69").Select 
Selection.Replace  What:="AB*",  Replacement:="AB6",  LookAt:=xlPait 

SearchOrder:=xlByRows,  MatchCase:=False 
Range("L9").  Select 

ActiveCell.FormulaRlCl  =  "Ln  TI*TiAB6" 
Range("L9").  Select 
Selection. Font.Bold  =  True 
With  Selection 

.  HorizontalAhgnment  =  xl  Center  AcrossSelection 

Vertical  Alignment  =  xlBottom 

WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 
Range("L9").  Select 

ActiveCell.FormulaRlCl  =  "Ln  Tl  *  TiAB6" 
With  ActiveCell.Characters(Start:=l,  Length:=13).Font 

Name  =  "Arial" 

FontStyle  =  "Bold" 

.Size  =  10 

Strikethrough  =  False 

Superscript  =  False 

Subscript  =  False 

OutlineFont  =  False 

Shadow  =  False 

Underline  =  xlNone 
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Colorlndex  =  xl  Automatic 
End  With 

Range("L10").SeIect 

ActiveCell.FormulaRlCl  =  "=RC[-3]*RC|-1]M 
Range("L10").Select 
Selection. AutoFill  Destinauon:=Range("L10:L69"),  Type: 

xlFillDefault 
Range("L10:L69").Select 
Range("M64").  Select 
ActiveWindow.LargeScroll  Down:=-3 
ActiveWindow.SmallScroll  Down:=12 
Active  Window.  ScrollRow  =  28 
Active  Window.  ScrollRow  =  1 
Active  Window.  ScrollRow  =  61 
Range("K70").  Select 
ActiveCell.FormulaRlCl  -  "SUM" 
Range("K70").Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
Selection.  Copy 
Range("L70").  Select 
ActiveSheet.Paste 
Range("K71").  Select 
Application.  CutCopyMode  =  False 
ActiveCell.FormulaRlCl  =  "=SUM(R[-61]C:R|-2]C)" 
Range("L71").Select 

AcuveCell.FormulaRlCl  =  "=SUM(R[-61]C:R[-2]C)" 
Range("L72").Select 
With  Toolbars(8) 

.Left  =  57 

.Top  =  106 
End  With 

ActiveWindow.LargeScroll  Down:=-4 
Range("M9").Select 
ActiveCell.FormulaRlCl  =  "Ti^" 
Range("M9").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("M10").Select 

ActiveCell.FormulaRlCl  =  "=RC[-5]AR[-4]C[4]" 
Range("M10").Select 
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Selection. AutoFill  Destination:=Range("M10:M69"),  Type:=_ 

xJFillDefault 
Range("M10:M69").Select 
Selection.Replace  What:="AQ*",  Replacement—"^",  LookAt:=xlPart, 

SearchOrder:=xlByRows,  MatchCase "False 
ActiveWindow.SmallScroll  Down:=8 
Active  Window.  ScrollRow  =  33 
Active  Window.  ScrollRow  =  43 
Active  Window.  ScrollRow  =  48 
Active  Window.  ScrollRow  =  51 
ActiveWindow.  ScrollRow  =  57 
Range("M70*').Select 
ActiveCell.FormulaRlCl  =  "SUM" 
Range("M70").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("M71").  Select 

ActiveCell.FormulaRlCl  =  "=SUM(R[-61]C:R[-2]Q" 
Range("M72").  Select 
ActiveCell.FormulaRlCl  =  "mean" 
Range("M73").  Select 
ActiveCell.FormulaRlCl  = "" 
Range("M72").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

VerticalAlignment  =  xlBottom 

.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("M73").  Select 
ActiveCell.FormulaRlCl  =  "=R[-2]C/60" 
Range("M74").Select 
Active  Window.  ScrollRow  =  1 
ActiveWindow.SmallScroll  ToRight:=2 
Range("N9").Select 

ActiveCell.FormulaRlCl  =  "(TiAQ6-mean)A2" 
Range("N9").Select 
With  Selection 

HorizontalAlignment  =  xlCenter  AcrossSelection 

VerticalAlignment  =  xlBottom 

WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 

Selection.Font.Bold  =  True 
ActiveWindow.SmallScroll  Down:=59 
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ActiveWindow.LargeScroll  Down:=-4 

Range("N10").Select 

ActiveCeU.FormulaRlCl  =  "=(RC[-11-R(63]C[-1])A2" 

Range("N10").Select 

Selection. AutoFill  Destination:=Range("N10:N69"),  Type:=  _ 

xlFillDefault 
Range("N10:N69").Select 
Selecuon.Replace  What:="-M*)A2",  Replacement:="-M73)A2",  LookAt 

:=xlPart  SearchOrder:=xlByRows,  MatchCase:=False 
Range("09").  Select 

ActiveCeU.FormulaRlCl  =  "(Ln  Ti)A2" 
Range("09").Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

VerticalAlignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("O10").Select 
ActiveCeU.FormulaRlCl  =  "=RC[-6]A2" 
Range("0 10").  Select 
Selection.AutoFill  Destination:=Range("O10:O69"),  Type:=  _ 

xlFillDefault 
Range("O10:O69").Select 
Range("M70").  Select 
Selection.  Copy 
Range("N70").Select 
ActiveSheet.Paste 
Range("O70").Select 
ActiveSheet.Paste 
Range("N71").Select 
Application.  CutCopy  Mode  =  False 
ActiveCeU.FormulaRlCl  =  "=SUM(R[-61]C:R[-2]Q" 
Range("071").  Select 

ActiveCeU.FormulaRlCl  =  "=SUM(R[-61]C:R[-2]Q" 
Range("072").  Select 
ActiveWindow.SmallScroll  ToRight:=2 
ActiveWindow.LargeScroll  Down:=-3 
Active  Window.  ScrollRow=  1 
ActiveWindow.LargeScroll  ToRight:=-l 
With  Toolbars(8) 

Left  =  543 

Top  =  37 
End  With 
Range("B4").Select 

ActiveWindow.SmallScroll  ToRight:=6 
ActiveWindow.SmallScroll  Down:=60 
ActiveWindow.LargeScroll  Down:=-4 
ActiveWindow.LargeScroll  ToRight:=- 1 
Range("B4").Select 
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ActiveCell.FormulaRlCl  =_ 

"=((6/(PI0A2))*(R[67]C[13]-((R[67]C[7]A2)/60))*(l/59))A(-0.5)" 
Range("B4").  Select 
Selection.  Copy 

Active  Window.  SmallScroll  ToRight:=10 
Range("Q4").  Select 
ActiveSheet.Paste 
Application.  CutCopyMode  =  False 
ActiveCell.FormulaRlCl  =_ 

"=((6/(PI0A2))*(R[67]C[-2]-((R[67]C[-8]A2)/60))*(l/59))A(-0.5)" 
Range("Q5").Select 

ActiveWindow.LargeScroll  ToRight:=-l 
Active  Window.  Small  Scroll  ToRight:=-l 
Range("C6").  Select 
ActiveCell.FormulaRlCl  =  _ 

"=(60/RC[-l])+R[65]C[6]-(60*R[65]C[9]/R[65]C[8])M 
Range("C7").Select 

Active  Window.  LargeScroll  ToRight  :=- 1 
Range("B9").  Select 

ActiveCell.FormulaRlCl  =  "=(60/R[62]C[9])A(l/R[-3]C)H 
Range("B  10").  Select 
Active  Window.  SmallScroll  ToRight:— 1 
Range("B  11").  Select 
ActiveCell.FormulaRlCl  =  "=1/R[-2]C" 
Range("B  12").  Select 
ActiveWindow.  SmallScroll  ToRight:=l  1 
Range("P6").  Select 

ActiveCell.FormulaRlCl  =  "=(R[65]C[-3])/((R[65]C[-2])A(0.5))" 
Range("Q9").  Select 

ActiveWindow.LargeScroll  ToRight:=-2 
ActiveWmdow.  SmallScroll  ToRight:=9 
Range("Q10").Select 
ActiveCell.FormulaRlCl  =  "theta" 
Range("Q10").Select 
Selection.Font.Bold  =  True 
With  Selection 

Horizontal  Alignment  =  xl  Center 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
Range("Q9").  Select 

ActiveCell.FormulaRlCl  =  "=1/(((1/60)*R[62]C[-4])A(1/R[-3)C))" 
Range("Qll").Select 
ActiveCell.FormulaRlCl  =  "=1/R[-2]C" 
Range("Q  12").  Select 
End  With 


End  Sub 
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APPENDIX  C.  MACRO  100  ENGINES 


MACRO  TO  100  ENGINES 

thesis  Macro 

Macro  recorded  by  FERNANDO  GUIMIL 


Sub  thesisO 
Sheets  Add 

Sheets("Sheet2").  Select 
Range("Al:C2").Select 

ActiveCeU.FormulaRlCl  =  "MAXIMUM  LIKELIHOOD  ESTIMATOR" 
Range("Al:C2").SeIect 
Range("A2")  Activate 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAhgnment  =  xlCenter 

.VerticalAhgnment  =  xlBottom 

.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
With  Selection 

HorizontalAhgnment  =  xlCenterAcrossSelection 

Vertical  Alignment  =  xlBottom 

WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 

RangefDl").  Select 

Active  Window.  SmallScroll  ToRight:=10 
With  Toolbars(8) 

.Left  =  241 

Top  =  228 
End  With 

Active Window.LargeScroll  ToRight:=-l 
Range("Pl:Rl").  Select 

ActiveCeU.FormulaRlCl  =  "MODIFIED  MOMENT  METHOD" 
Range("Pl:Rl").  Select 
Range("Ql").  Activate 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAhgnment  =  xlCenter 

VerticalAhgnment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
With  Selection 

.HorizontalAhgnment  =  xlCenterAcrossSelection 

VerticalAhgnment  =  xlBottom 
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.  WrapText  =  True 

.Orientation  =  xlHorizontaJ 
End  With 
Range("Nl").Select 

Active Wi ndow. LargeScroll  ToRight :  =- 1 
With  Toolbars(8) 

Left  =  436 

Top  =  226 
End  With 
Range("A2").Select 

ActiveCell.FormulaRlCl  =  "shape  parameter" 
With  ActiveCell.Characters(Start:=l,  Length:=15).Font 

.Name  =  "AriaJ" 

FontStyle  =  "Bold" 

.Size  =10 

Strikethrough  =  False 

Superscript  =  False 

Subscript  =  False 

OutlineFont  =  False 

Shadow  =  False 

.Underline  =  xlNone 

Colorlndex  =  xl  Automatic 
End  With 
Range("B3").  Select 

ActiveCell.FormulaRlCl  =  "initial  beta" 
Range("B3").  Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

Vertical  Ah  gnment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
Range("B5").Select 

ActiveCell.FormulaRlCl  =  "final  beta" 
Range("B5").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("A7:C7").  Select 
ActiveCell.FormulaRlCl  =  "scale  parameter" 
Range("A7:C7").Select 
Range("B7").Activatc 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 
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.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("A7").  Select 
ActiveCell.FormulaRlCl  =  "" 
Range("B7").  Select 

ActiveCell.FormulaRlCl  =  "scale  parameter" 
Range("B8").  Select 
ActiveCell.FormulaRlCl  =  "landa" 
Range("B8").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

AcuveWindow.  Small  Scroll  ToRight:=0 
Range("A2:C2").Select 
Selection.  Copy 

Active  Window.  SmallScroll  ToRight:=l  1 
Range("Q2").Select 
ActiveSheet.Paste 
Application.  CutCopy  Mode  =  False 
Selection.  Clear 

ActiveWindow.LargeScroll  ToRight:=-2 
Range("A2:C2").Select 
Selection.  Copy 

AcuveWindow.  SmallScroll  ToRight:=9 
Range("P2").  Select 
ActiveSheet.Paste 
Range("04").Select 

ActiveWindow.LargeScroll  ToRight:=- 1 
Range("B3").  Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 

Active  Window.  SmallScroll  ToRight:=l  1 
Range("Q3").  Select 
ActiveSheet.Paste 

ActiveWindow.LargeScroll  ToRight:=-2 
Range("B5").  Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 

Active  Window.  SmallScroll  ToRight:= 1 2 
Range("Q5").Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 
ActiveSheet.Paste 
Range("R7").  Select 

Active  Window.  LargeScroll  ToRight : =-2 
Range("B5").Select 
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Application.  CutCopyMode  =  False 

Selection.  Copy 

ActiveWindow.  Small  Scroll  ToRight:=12 

Range("Q5").Select 

ActiveSheet.Paste 

Active  Window.LargeScrollToRight:=-2 

Range("A7:C8").Select 

Application.CutCopyMode  =  False 

Selection.  Copy 

ActiveWindow.  SmallScroll  ToRight:=l  1 

Range("P7").  Select 

ActiveSheet.Paste 

Range("N7").Select 

ActiveWindow.LargeScroll  ToRight:=-2 

Range("D  11").  Select 

Application.CutCopyMode  =  False 

Selection.  Copy 

Range("D7").  Select 

ActiveSheet.Paste 

Range("B8").  Select 

ActiveWindow.  Small  Scroll  ToRight:=13 

ActiveWindow.LargeScroll  ToRight:=-2 

ActiveWindow.  SmallScroll  ToRjght:=14 

ActiveWindow.LargeScroll  ToRight:=-2 

Range("El").  Select 

Application.CutCopyMode  =  False 

ActiveCell.FormulaRlCl  =  "Number  of  observations" 

Range("El").Select 

Selection.Font.Bold  =  True 

With  Selection 

Horizontal  Alignment  =  xl  Center 

VerticalAlignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
With  Selection 

HorizontalAlignment  =  xlCenterAcrossSelection 

VerticalAlignment  =  xlBottom 

WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 
Range("E2").Select 
ActiveCell.FormulaRlCl  =  "n" 
Range("E2").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

VerticalAlignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
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With  Toolbars(8) 

Left  =  539 

Top  =  30 
End  With 

Sheets("Sheetl").Select 
Range("B2:B  101").  Select 
Selection.  Copy 
Sheets("Sheet2").  Select 
Acuve  Window.  SmallScroll  Down:=l 
Range("D  10").  Select 
AcuveSheet.Paste 
SheetsfSheetl").  Select 
Range("C2:C101").Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 
Sheets("Sheet2").  Select 
Range("E10").Select 
AcuveSheet.Paste 
Sheets("Sheetl").  Select 
Range("D2:D101").Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 
Sheets("Sheet2").Select 
Range("F10").Select 
ActiveSheet.Paste 

AcuveWindow.LargeScroll  ToRight:=0 
Range("D10:D109").Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 
Range("G10").Select 
AcuveSheet.Paste 
Range("E10:E109").Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.SmallScroll  Down:=100 
Range("G110").Select 
ActiveSheet.Paste 

Acuve  Window.  LargeScroll  Down:=-6 
Range("F10:F109").Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 

AcuveWindow.SmallScroll  Down:=200 
Range("G210").Select 
ActiveSheet.Paste 

Acuve  Window.  LargeScroll  Down:=-l  1 
Range("G10:G309").Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 
Range("H10").Select 
AcuveSheet.Paste 
ApplicaUon.  CutCopyMode  =  False 
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Selection.  Sort  Keyl:=Range("H10"),  Order l:=xl Ascending,  Header: 

xlGuess,  OrderCustom:=l,  MatchCase-False,  Orientation— _ 

xlTopToBottom 
ActiveWindow.SmallScroll  Down:=-l 
Range("H2:H6").Select 
ActiveCell.FormulaRlCl  =  "failures  times  Ti" 
Range("H2:H6").Select 
Range("H3").  Activate 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
With  Selection 

.HorizontalAlignment  =  xlCenter  AcrossSelection 

.VerticalAlignment  =  xlBottom 

.WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 

Range("H9").  Select 
ActiveCell.FormulaRlCl  =  "Ti" 
Range("H9").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

.VerticalAlignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
Range("I9").  Select 
ActiveCell.FormulaRlCl  =  "Ln  Ti" 
Range("I9").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

VerticalAlignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Rangefl  10").  Select 

ActiveCell.FormulaRlCl  =  "=LN(RC[-1])M 
Range("I10").Select 
Selection.  AutoFill  Destination:=Range("1 10:1309"),  Type:=  _ 

xlFillDefault 
Range("110:I309").Select 
Willi  Toolbars(8) 

Left  =  79 

Top  =  270 
End  With 
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Range("G315").Select 

Active  Window.  SmallScroll  ToRight:=2 

Range("I305").Select 

ActiveWindow  .LargeScroll  Do\vn:=- 1 7 

Range("J9").  Select 

Active  Window. LargeScroll  ToRight:=-l 

ActiveWindow.  Small  Scroll  ToRight:-2 

Range("J9").  Select 

ActiveCell.ForrnulaRlCl  =  ,,TiAB4" 

Range("J9").  Select 

Selection.Font.Bold  =  True 

With  Selection 

.HorizontalAlignment  =  xJCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
Selection.  Copy 
Range("K9").  Select 
ActiveSheet.Paste 

ActiveWindow. LargeScroll  ToRight : =- 1 
ActiveWindow. SmallScroll  ToRight:=4 
Application.  CutCopyMode  =  False 
ActiveCell.ForrnulaRlCl  -  "TiAB6" 
Range("J10").Select 

AcuveCeU.FormulaRlCl  =  "=(RC[-2]AR[-6]C[-8])M 
Range("J  10").  Select 
Selection. AutoFill  Destination:=Range("J10:J309"),  Type:=  _ 

xlFillDefault 
Range("J10:J309").Select 
Range("J300").  Select 
ActiveWindow.LargeScroll  Down:=-16 
ActiveWindow.SmallScroll  Down:=-9 
Range("J10:J309").Select 
Selection.Replace  What:- 'AB*)'\  Replacement:=::"AB4)",  LookAt: 

xlPart,  SearchOrder:=xlByRows,  MatchCase:=False 
Range("K10").Select 

ActiveCell.ForrnulaRlCl  =  M=(RC[-3]AR[-4]C[-9])" 
Range("K  10").  Select 
Selection.  AutoFill  Destination:=Range("K10:K309"),  Type:=_ 

xlFillDefault 
Range("K10:K309").Select 
SelecUon.Replace  What:="AB*)",  Replacement:="AB6)",  LookAt: 

xlPart,  SearchOrder:=xlByRows,  MatchCase "False 
Range("L8:L9").Select 
ActiveCell.FormulaRlCl  =  "Ln  Ti  *  TiAB6" 
Range("L8:L9").Select 
Range("L9").Activate 
With  Selection 

HorizontalAlignment  =  xlCenterAcrossSelection 

VerticalAlignment  =  xlBottom 
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WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 
Range("L9").Select 
Selection.Delete  Shift:=xlUp 
Range("L8").  Select 
Selection.  Copy 
Range("L9").  Select 
ActiveSheet.Paste 
Range("L8").  Select 
Application.  CutCopyMode  =  False 
Selection.  Clear 
Range("L9").  Select 
Selection.  Font.  Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

.VeiticalAlignment  =  xJBottom 

WrapText  -  True 

Orientation  =  xlHorizontal 
End  With 

Range("L  10").  Select 

ActiveCell.FormuIaRlCl  =  "=RC[-3]*RC[-1]" 
Range("L10").Select 
Selection. AutoFill  Destination:=Range("L10:L309"),  Type: 

xlFillDefault 
Range("L10:L309").Select 
ActiveWindow.LargeScroll  Down:=-l 
Range("M285").  Select 
AcuveWindow.LargeScroll  Down:=-16 
Active  Window.  SmallScroll  ToRight:=4 
ActiveWindow.LargeScroll  ToRight:=-l 
ActiveWindow.  SmallScroll  ToRight:=8 
Active  Window.  SmallScroll  Down:=33 
ActiveWindow.LargeScroll  Down: =-2 
ActiveWindow.  ScrollRow  =  174 
ActiveWindow.  ScrollRow  =  218 
ActiveWindow.  ScrollRow  =  252 
ActiveWindow.  ScrollRow  =  309 
Range("L3 10").  Select 
ActiveCell.FormuIaRlCl  -  "SUM" 
Range("L3 10").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("L3 11").  Select 

ActiveCell.FormuIaRlCl  =  "=SUM(R[-301]C:R[-2]Q" 
Range("L312").Select 
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Range("M9").  Select 
ActiveCell.FormulaRlCl  =  "Ti^" 
Range("M9").  Select 
Selection.Font.Bold  =  True 
With  Selection 

.Horizontal  Alignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("M10").Select 

ActiveCell.FormulaRlCl  -  "=(RC[-5]AR[-4]C[4])" 
Range("M10").Select 
Selection. AutoFill  Destination:=Range("M10:M309").  Type:=  _ 

xlFillDefault 
Range("M10:M309").Select 
Selection.Replace  What:=,,AQ*)",  Replacement:^"^)".  LookAt: 

xlPart,  SearchOrder:=xlByRows,  MatchCase:=False 
Range("N9").  Select 

Acuve Window. SmallScroll  Down:=107 
Active  Window.  ScrollRow=  162 
Active  Window.  ScrollRow  =  214 
Active  Window.  ScrollRow  =  243 
ActiveWindow.  ScrollRow  =  261 
Active  Window.  ScrollRow  =  282 
ActiveWindow.  ScrollRow  =  297 
Range("M310").Select 
AcuveCell.FormulaRlCl  =  "SUM" 
Range("M310").Select 
Selecuon.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

VerticaLAlignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("M311").Select 

AcuveCell.FormulaRlCl  =  "=SUM(R[-301]C:R[-2]Q" 
Range("M312").Select 
ActiveCell.FormulaRlCl  =  "mean" 
Range("M312").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

.Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("M3 13").  Select 
ActiveWindow.LargeScroll  Down:=-17 
Range("H3").  Select 
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Active Window.LargeScroll  ToRight  :=- 1 
Range("E3").  Select 
ActiveCell.FormulaRlCl  =  "300" 
Range("E3").  Select 
Selection. Font. Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Active  Window.  SmallScroll  ToRight:=5 
ActiveWindow.ScrollRow  =  144 
Active  Window.  ScrollRow  =  178 
ActiveWindow.ScrollRow  =  219 
Active  Window.  ScrollRow  =  264 
ActiveWindow.ScrollRow  =  291 
ActiveWindow.ScrollRow  =  312 
ActiveWindow.ScrollRow  =  274 
ActiveWindow.LargeScroll  Down:=-15 
ActiveWindow.  SmallScroll  Down:=37 
ActiveWindow.ScrollRow  =  283 
ActiveWindow.ScrollRow  =  312 
ActiveWindow.ScrollRow  =  303 
Range("M313").Select 
ActiveCell.FormulaRlCl  =  "=R[-2]C/300" 
Range("M314").Select 
ActiveWindow.LargeScroll  Down:=-17 
Range("N9").  Select 

ActiveCell.FormulaRlCl  =  "(TiAQ6-mean)A2" 
Range("N9").  Select 
Selection. Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

VerticalAlignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
With  Selection 

HorizontalAlignment  =  xlCenterAcrossSelection 

VerticalAlignment  =  xlBottom 

WrapText  =  True 

Orientation  =  xlHorizontal 
End  With 

Range("N  10").  Select 

ActiveCell.FormulaRlCl  =  "=(RC(-2]-R[303]C[-l])H 
Range("N10").Select 
Selection. AutoFill  Destination:=Range("N10:N309").  T>pe: 

xlFillDefault 
Range("N10:N309").Select 
ActiveWindow.SmallScrollToRight:=l 
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Range("N10").Select 

ActiveCell.FormulaRlCl  -  "=(RC[-1]-R[303|C[-1|)A2" 

Range("N10").Select 

Selection.  AutoFill  DestinaUon:=Range("N10:N309"),  Type:=  _ 

xlFillDefault 
Range("N10:N309").Select 
Selection. Replace  What:="-M*)A2",  Replacement:="-M313)A2",  LookAt 

:=xlPart  SearchOrder:=xlByRows.  MatchCase:=False 
ActiveWindow.SmallScroll  Down:=-2 
ActiveWindow.SmallScroll  ToRight:=3 
ActiveWindow.LargeScroll  ToRight:=-l 
Range("B4").  Select 

ActiveCell.FormulaRlCl  =  "=((6/(PI0A2))*l)" 
Range("B5").Select 

ActiveWindow.SmallScroll  Down:=16 
Active  Window.  ScrollRow=  1 
ActiveWindow.SmallScroll  Down:=20 
Active  Window.  ScrollRow  =  612 
Active  Window.  ScrollRow  =  581 
Active  Window.  ScrollRow  =  539 
ActiveWindow.  ScrollRow  =  452 
ActiveWindow.LargeScroll  Down:=-9 
ActiveWindow.  ScrollRow  =  339 
ActiveWindow.  ScrollRow  =  309 
Range("H3 10").  Select 
ActiveCell.FormulaRlCl  =  "SUM" 
Range("H3 10").  Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("H3 11").  Select 

AcuveCell.FormulaRlCl  =  "=SUM(R[-301]C:R[-2]Q" 
Range("H3 12").  Select 
ActiveWindow.  ScrollRow  =  1 
ActiveWindow.  ScrollRow  =  465 
ActiveWindow.  ScrollRow  =  287 
ActiveWindow.  ScrollRow  =  322 
ActiveWindow.LargeScroll  ToRight:=l 
Active  Window.  ScrollRow  =  203 
ActiveWindow.  ScrollRow  =  269 
ActiveWindow.  ScrollRow  =  3 
ActiveWindow.SmallScroll  ToRight:=-4 
Active  Window.  ScrollRow  =  124 
ActiveWindow.  ScrollRow  =  223 
ActiveWindow.  ScrollRow  =  317 
ActiveWindow.  ScrollRow  =  274 
ActiveWindow.  ScrollRow  =  276 
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Active  Window.  SmallScroll  Down:=23 

Range("H310").Select 

Selection.  Copy 

Range("I310").Select 

ActiveSheet.Paste 

Range("I3 11").  Select 

Application.  CutCopy  Mode  =  False 

ActiveCell.FormulaRlCl  =  "=SUM(R[-301]C:R[-2]C)" 

Range("I312").Select 

Active  Window.  ScrollRow=  1 

Active  Window.  SmallScroll  Down:=l 

Active  Window.  SmallScroll  ToRight:=2 

Range("09").  Select 

ActiveCell.FormulaRlCl  =  "(Ln  Ti)A2" 

Range("09").Select 

Selection.Font.Bold  =  True 

With  Selection 

.HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("O10").Select 
ActiveCeU.FormulaRlCl  =  "=RC[-5]A2" 
Range("O10").Select 
Selection.  AutoFill  Destination:=Range("O10:O309"),  Type:=  _ 

xlFillDefault 
Range("O10:O309").Select 
Range("03 10").  Select 
ActiveCell.FormulaRlCl  =  "SUM" 
Range("03 10").  Select 
Selection.Font.Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  Willi 

Range("03 11").  Select 

ActiveCell.FormulaRlCl  =  "=SUM(R[-301]C:R[-2]Q" 
Range("03 12").  Select 
ActiveWindow.  Small  Scroll  ToRight:=-7 
Active Window.LargeScroll  Down:— 15 
Active  Window.  ScroIlRow=  1 
Active  Window.  Small  Scroll  ToRight:=l 
Range("B4").  Select 
ActiveCell.FormulaRlCl  =_ 

"=((6/(PI0A2))*(R[307]C[13]-R(307]C[7]/300)/299)A(-0.5)" 
Range("B4").  Select 

ActiveWindow.  SmallScroll  ToRight:=6 
ActiveWindow. SmallScroll  Down:=-l 
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ActiveCell.FormulaRlCl  =  _ 

"=((6/(PIOA2))*(R[307]C|13)-(R|307]C[7]A2)/300)/299)A(-0.5)" 
Range("B4").Select 
Active  Window.  ScrollRow  =  419 
Active  Window.  ScrollRow  =  300 
ActiveWindow.SmallScrollToRight:=7 
ActiveWindow.ScrollRow  =  1 
Active Window.LargeScroll  ToRight.-- 1 
ActiveCell.FormulaRlCl  =_ 

"=((6/(PI0A2))*(R[307]C[13]-(R[307]C[7]A2)/300)/299)A(-O.5)M 
Range("B4").  Select 
ActiveCell.FormulaRlCl  =_ 

"=((6/(PI0A2))*(R[307]C[13]-(R[307]C[7]A2)/300)/299)A(-0.5)" 
Range("1 10").  Select 
ActiveWindow.  Small  Scroll  ToRight:=2 
Range("O10").Select 
ActiveWindow.  LargeScroll  ToRight:=-l 
Range("B4").  Select 
ActiveCell.FormulaRlCl  =_ 

"=((6/(PI0A2))*(R[307]C[13)-(R[307]C[7]A2)/3OO)/299)A(-0.5)M 
Range("B4").Select 

ActiveWindow. SmallScroll  ToRight:=2 
Range("J10").Select 

ActiveWindow.  LargeScroll  ToRight  :=-l 
ActiveWindow. SmallScroll  ToRight:=2 
ActiveWindow.LargeScroll  ToRight:=-l 
Range("B4").  Select 

ActiveWindow.SmallScroll  ToRight:=2 
Range("Jll").Select 

ActiveWindow.LargeScroll  ToRight:=- 1 
Range("B4M).  Select 
ActiveCell.FormulaRlCl  =_ 

"=((6/(PI0A2))*(R[307]C[13]-(R[307]C[7]A2)/300)/299)A(-0.5)" 
Range("O10").Select 
ActiveWindow.LargeScroll  ToRight:=-l 
ActiveWindow.  SmallScroll  ToRight:=7 
ActiveCell.FormulaRlCl  =  "=RC[-6]A2" 
Range("O10").Select 
Selection.AutoFill  Desunation:=Range("O10:O309''),  Type:=_ 

xlFillDefault 
Range("O10:O309").Select 
ActiveWindow.LargeScroll  Down:=-16 
AcuveWindow.SmallScroll  Down:=-12 
ActiveWindow.LargeScroll  ToRight:=-l 
Range("B4").Select 

Active  Window.  SmallScroll  ToRight:=3 
ActiveWindow.LargeScroll  ToRight:=-l 
AcuveWindow.SmallScroll  ToRight:=13 
Range("Q4").  Select 

ActiveCell.FormulaRlCl  =  "=RC[-15]" 
Range("Q5").Select 
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ActiveWindow.LargeScrollToRight:=-2 
With  Toolbars(8) 

Left  =  398 

Top  =139 
End  With 

ActiveWindow.LargeScroll  ToRight:= 1 
ActiveWindow.SmallScroll  ToRight:=-9 
Active Window. SmallScroll  Down:=3 
Range("C6").Select 
ActiveCell.FormulaRlCl  =_ 

"=(300/RC[-l])+R[305]C[6]-(300*R[305]C[9]/R[305]C[8])" 
Range("K310").Select 
ActiveCell.FormulaRlCl  =  "SUM" 
Range("K310").Select 
Selection. Font. Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("K311").Select 

ActiveCell.FormulaRlCl  =  "=SUM(R[-301]C:R[-2]C)" 
Range("K312").Select 
ActiveWindow.LargeScroll  Down:=-17 
ActiveWindow.LargeScroll  ToRight:=- 1 
Range("C6").  Select 
ActiveCell.FormulaRlCl  =_ 

"=(300/RC[-l])+R[305]C[6]-(300*R[305]C[9]/R[3051C[8])" 
Range("C6").  Select 

ActiveWindow.SmallScroll  ToRight:=2 
Active  Window.  ScrollRow  =  432 
Active  Window.  ScrollRow  =  311 
Active  Window.  ScrollRow  =  1 
Active  Window.  S  mall  Scroll  ToRight : =-2 
Range("B10").Select 
ActiveCell.FormulaRlCl  =  "theta" 
Range("B  10").  Select 
Selection.  Font.  Bold  =  True 
With  Selection 

HorizontalAlignment  =  xlCenter 

Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 
Range("B9").  Select 

ActiveCell.FormulaRlCl  =  "=(300/R[3021C[9])A(l/R[-3]C)" 
Rangc("B9").  Select 
ActiveCell.FormulaRlCl  = "" 
Range("C9").  Select 
ActiveCell.FormulaRlCl  =  "" 
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Range("B9").  Select 

ActiveCell.FormulaRlCl  =  "=(300/R|302|C|9])A(l/R[-3]C)" 

RangeC'B  11").  Select 

ActiveCell.FormulaRlCl  =  "=1/R[-2]C" 

Range("B  II").  Select 

Active  Window.  SmallScroll  ToRight:=9 

Range("Q6").  Select 

With  Toolbars(8) 

Left  =  58 

Top  =  25 
End  With 
Range("P6").Select 

ActiveCell.FormulaRlCl  =  "((1/299)" 
Range("S319").Select 
Active  Window.  ScrollRow  =  302 
Range("N310").Select 
ActiveCell.FormulaRlCl  =  "SUM" 
Range("N310").Select 
Selection.Font.Bold  =  True 
With  Selection 

.  Horizontal  Alignment  =  xl  Center 

.Vertical  Alignment  =  xlBottom 

.WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

Range("N311").Select 

ActiveCell.FormulaRlCl  -  "=SUM(R[-301]C:R[-2]C)M 
Range("N312").Select 
Active  Window.  LargeScroll  Down:=-17 
Range("P6").  Select 

ActiveCell.FormulaRlCl  =  "(M311)/((N311)A(0.5))" 
Range("P6").  Select 

Active  Window.  LargeScroll  ToRight:=0 
ActiveCell.FormulaRlCl  =_ 

"=R[305]C[-3])/(R[305]C[-2])A(0.5)" 
Range("Q9").Select 

ActiveCell.FormulaRlCl  =  "=1/(((1/300)*R[302]C[-41)A(1/R[-3]C))" 
Range("Q  10").  Select 
Active Window.LargeScroll  ToRight:=- 1 
AcUveWindow.  Small  Scroll  ToRight:=13 
ActiveCell.FormulaRlCl  =  "theta" 
Range("Q10").Select 
Selecuon.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 

.Vertical  Alignment  =  xlBottom 

WrapText  =  False 

Orientation  =  xlHorizontal 
End  With 

RangeC'Q  11").  Select 
ActiveCell.FormulaRlCl  =  "=1/R[-2]C" 
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Range("Q12').Select 
ActiveWindow.LargeScroll  ToRight:=-2 
End  Sub 
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APPENDIX  D.  PARAMETER  CALCULATIONS  FOR  THE  20  AND  100  ENGINES 


Row  Names 
1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 


V1 

V2 

V3 

88.980852932 

7.374101636 

166.612433434 

57.770945753 

452.164077047 

671.533220087 

348.363054836 

13.725217444 

33.719418711 

382.257870284 

168.795667771 

119.352260671 

57.145359689 

117.795729216 

228.544288266 

251.488400993 

417.099715694 

206.169150245 

703.336492627 

16.866142294 

168.320579106 

85.174113550 

960.378668893 

561.250692879 

126.379607253 

166.962891751 

410.881188239 

9.871149363 

30.987527462 

1171.995737937 

184.016586935 

55.183363300 

170.715797827 

19.085994032 

206.177991742 

124.656084035 

171.708369596 

7.411273439 

458.366290751 

957.877021060 

334.941605556 

252.553811069 

88.483061291 

129.532846338 

55.052226186 

97.662694658 

64.764644477 

521.223680534 

213.546020387 

8.955316669 

314.364020817 

420.340591055 

182.036946622 

345.506361499 

46.823610226 

411.857588415 

109.450083602 

25.300667331 

361.149199434 

744.757337329 

Column  1            Column  2            Column  3 

Column  1 
Column  2 
Column  3 

1 

-0.099201099 

1 

-0.325682845 

0.177625785 

1 
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ABC 

D 

E 

f 

H 

1                   J 

K 

L 

M 

N 

O 

P                Q                R 

1 

MAXIMUM  LIKELIHOOD 
ESTIMATOR 

Number 

ot 

observati 

ons 

Failures 
times  Ti 

MODIFIED  MOMENT  METHOD 

2 

shape  parameter 

n 

shape  parameter 

3 

initial  beta 

60 

initial  beta 

4 

1  001851 

1001851 

5 

tinal  beta 

final  beta 

6 

0973 

-0  00065 

7811979 

0  985725 

7 

scale  parameter 

scale  parameter 

8 

land.! 

landa 

9 

0  003973 

Ti 

LnTi 

TiAB4 

|    LnTV 
TiAB6         Ti"B6 

Ti'QS 

(Trot 

mean) '2 

<Ln  Ti)*2 

000395 

10 

theta 

88  98085 

7374102 

1666124 

88  98085 

7  374102 

1.997974 

740142 

6.986841,    13.95953!  7.166757 

5141343        3.9919 

theta 

11 

251  7053 

57  77095 

452  1641 

671  5332 

57  77095'   7  411273 

2  003002 

7  438798 

8  991725 

9  913067 

70211081    1406329,   7  202367    5139729 
8  440626J    18  50394    8679408    50729  75 

4  012018 

2531397 

12 

348  3631 

13  72522     33  71942     348  3631 
"  168  7957*"  119  3523'   382  2579 

8  955317,   2  192247 

4  805949 

13 

3822579 

9  8711491  2.289616 

9  279395    2124625 

9.553737    50336  66 

5  242343 

14 

57  14536 

117  7957 

228  5443    57  14536 
206  1692    251  4884 

13  72522,   2619235 

1379191 

12.7881 

3349504 

13.22152 

48704  32 

47398  78 

"46488  89 

43995  8 

41782  11 
40741  09 

6  860391 

15 

251  4884 

417  0997 

16  86614    2  825308 

1695456 

1562738 

4415217 

16  19945 

18  29923 

7  982366 

16 

703  3365 

16  86614     168  3206 

703  3365 

19  08599^2.948955 

19  19045 

17  62525 

51  97607 

8696334 

17 

85174111  9603787'  5615507 

85  17411 

25  30067I   3230831 

25  4524 

23  18714 

74  91374 

24  16031 
29  50533 

10  43827 

18 

126  3796     1669629    4108812 

1263796 

30  98753    3.433585 

31  18507 

28  2439 
30  66388 

96  97783 

11  7895 

19 

9871149    3098753    1171996 

9  871149 

33  71942,   3518074    33  93968 

107  8778    32  06785 

12  37684 

20 

184  0166,   5518336    170  7158    184,0166 

46  82361!   3  8463881  4715812 

42  20483    162.3361    44  32199      35944  4 
49  40532*  198  0305    5199069    33095  39 
49  51 983  "  1 98  6073    52  1 1 276*  33050  99 

147947 

21 

19  085991     206.1781   124.6561     19.08599J  55  05223 

^4  008282^55  46214 
*~4  010662      55  5945 
*T  045598    57  57483 

16  06633 

22 

171  7084 

7411273 
334  9416 

458  3663    1717084155  18336 

16  08541 

23 

957  877 

252  5538      957  877    57  14536 

51  23211 
51  77774 

207  2645 

53  93866 

32390  42 
32181  27 
29891  22 

16  36686 

24 

88  48306    129  5328 

55  05223 

88  48306 

97  66269 

57  77095    4  056486 

58  20629 
65  2665 

210,0357 

54  52069 

16  45508 

25 

97  66269 

64  76464    521 .2237 

64.76464|     4  17076 

57  86709 

241.3497 

61.02129 

17  39524 

26 

213  546 

8955317      314364      213  546    85  1741 1 
182  0369*  345  5064*"420  34061  88  48306 

4  444698    85  87764 
4  482811    89  22022 

755421 

78  39614 
78  82524 

86  299 

335  7618      79.9379 

23708  04 
22774  98 

19  75534 

27 

420  3406 

3514351     82  99826 

20  0956 

28 

46  82361 

411  8576*   109  4501 

46  82361!   88  98085 

4  488421 

89  72309 

98  49432 
110  4054 

353  8009 

83  45851 

22636  27 

20  14592 

29 

25  30067 

361  1492 

744  7573 

25  30067!  97.66269 

4  58152 

395  3806 

91  47988^  20286  93 
102  35441   17307  42 

20  99032 

30 

7  374102 

109  4501 

4  695469 

96  41775|  452.7265 

22  04743 

31 

452  1641 

1177957J  4768952|       11884 

1193523    4  782079'    120  4133 

""124  6561,   4  8255591   125.7744 

103  564 

104  8953 

493  8917 

110  0435    15343  42      22  7429 

32 

13  72522 

5016176,   1114767    1499041 

22  86828 

33 

168  7957 

109  4281     528  0518     116  3583     13818  88 

23  28602 

34 

117  7957 

1263796,     483929,   1275166 

110  9 
113  5914 

536  6772      117  944     13448  59 

2341873 

35 

417  0997 

1295328 

4  863934 

130  7042 

552  5011'   120  8442    12784  33!  23  65786 

36 

1686614 

1666124 

5.11567 

1681974 

145  11791   742  3754    154  8791    6246  208 

26  17008 

37 

9603787    1669629 

5.117772 

1685518 

14541491    7442003    1552002 
146  5653    7512749    156  4442 

6195552 

26  19159 

38 

166  9629    168  3206      5  12587 

169  925 

6001  273 

26  27455 

39 

30  98753     168  7957     5  128689 

170  4055 

146  9678 

753  7523     156  8794 

5934  025 

26  30345 

40 

55  18336     170  7158             5  14 

172  3475 

148  5943 

763  7745     158  6384     5666  125 

264196 

41 

206  178    1717084    5  145798:    173  3515 

1494348    768  9614    159  5475 

5530  082 
4213  102 

3981  217 
183"5  595 

26  47923 

42 

7  411273    182  0369 
334  9416'    184  0166 

5  20421;    183  7987, 

158  1739j  823  1704 

169  0036 

170  8152 

27  0838 

43 

5215026'   1858013 

159  8474'    8336083 

27  1965 

44 

129  5328 

206  1692|    5  328697J   208  2125 
206  178       5  32874,   2082214 

178  5416    951  3939,   191  0682 

28  39501 

45 

64  76464 

178549    9514413    1910763    1834  903 
184  7545     990  9956     197  8055*   1303  685 
197  3685*~1072O53~211  4932    502  6057 

28  39547 

46 

8  955317 

213546J  5363852 

215  6765 
*"  230  8534 

28  77091 

47 

182  0369 

228  5443,     5  43173 

29  50369 

48 

411  8576 

251  4884    5  527397    254  0743 

252  5538    5  531624    2551526 

'     314  364    5  750552    317  7276 
334  9416    5.813956    338  5651 
345  5064^5  845011     349  2642 
348  3631*  5  853245~T52  1573 
361  1492i  5  8892911  3651071 
382  2579'    5  946095'    386  4877 

216  6226 

217  5154 

1197  359    232  4079 

2  262488 
0  284797 
3099  941 

30  55212 

49 

361  1492 

1203  214    233  3784 

30  59887 

50 

1666124 

269  1546     1547  788    289  5892 

33  06884 

51 

671  5332 

286  2824     1664  433    308  2659 
295  0648'    1724~'657*"3178483 
297  4383'    1740  979    320  4387 
308  0554'   1814  228'     332  029 
325.5611     1935.817    351  1508 

349  25741    2101  937     "377  056 

350  06491  2107  626*3779392 
354  3995J  2138  207*  382  6805 

5528  496 
7045  298 
7486  855 
9626  936 
13744  92 
20490  18 
20743  81 
22132  05 
2301267 

33  80209 

52 

33  71942 

34  16415 

53 

119  3523 

34  26048 

54 

228  5443 

34  68375 

55 

206  1692 

35  35605 

56 

168  3206 

410  8812    6  018304    415  4833 

36  21998 

57 

561  2507 

411  8576,   6  020678 

416  4725 
421  7832 

36  24856 

58 

4108812 

417  0997'  6  033325 

36  40101 

59 

1 1 7 1  996 

420  34061  6  041065    425  0666 

357  0786 
383  3564 
388  4719 

2157  135    3856113 

36  49447 

60 

1707158 

452  1641     6  114045    457  3096 
'  458  3663^  6  127669-  463  5941 

2343  858    414  3736 

32566  36    37  38155 

61 

124  6561 

2380  427 

4199758 

476  6932 

"    512  758" 

34619  7     37  54832 

63 

252  5538 

5612507'  6  330168'   567  8647 
*  671  5332*  6  509563'  679  6725 

473  0743'    2994  639 

'   77755  46 

40  07102 

64 

5505223 

563  2956    3666  808    6119437       142908 

42  37442 

65 

5212237     703  3365    6  555835    7119222 

5892362     3862  935    6405017     165315  1     42  97898 

66 


A 

B                  C 

D 

E 

F                     G                     H                      1 

J                  K                  L                  M 

NOP 

Q 

66 

314  3641   744  75731  6  613058 

753  9285I  622  9742 

4119  7651  677  6683 

196919  61    43  73254 

67 

345  5064'     957  877    6  864719'   970.12431   7958185 

5463  071     868.4643 

402656.5!  47  12437 

68 

109  4501 

960  3787 

6  867328    972  6626 

797  8407 

5479.034           870.7 

4054989'   47  16019 

69 

744  7573 

1171  996 

7066463 

1187  424 

968  422 

68433191     105954 

6816622 

49  9349 

70 

SUM 

SUM 

SUM 

SUM            SUM 

SUM 

SUM 

71 

15284  6 

296  1914 

13008  26 

77584.93     14034.72 

3227644 

1558  848 

7  7 

mean 

73 

| 

233.9121| 

67 


Row  Names 
1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 


V1 

V2 

V3 

450  634670353 

134  620284266 

763.602013650 

395264314931 

300696009190 

201  433889226 

467  906036969 

301  510351673 

22  166715862 

200004796884 

261  419709923 

234019837022 

62  187250449 

351.566316602 

293954533661 

188420555520 

530.441439163 

446  499347177 

517  045829232 

212.650323103 

4  372290612 

104  406621033 

91  440692508 

286  231002289 

185580658025 

189  658931128 

461  452317203 

278  805459607 

142  687073420 

63.785901768 

242.374076871 

140  003133652 

57  571537504 

191  638384699 

35  982496826 

135685578893 

115.651429115 

807515319374 

56.217725189 

318.738490247 

220.511514630 

96  958868093 

389.324028552 

121  699393963 

391  967785281 

108.549643630 

833  642731237 

587  010089214 

647  283140422 

270  006356536 

600.593877204 

135.277879814 

235  260576261 

294.610909130 

167  565552866 

19.382052477 

80.687152199 

620417102151 

43940840941 

348.339229735 

549652520297 

147263663645 

402.329146256 

160858623029 

198697751449 

131  936073250 

269748636292 

216407382425 

272564254146 

119  882883842 

42753769287 

568  105938757 

290  367131217 

46018913125 

470  135356798 

111  526180647 

59.040162467 

442.572648090 

315.869668084 

112.068427109 

799.073849390 

258473619251 

740  012832046 

730  391096774 

936074568617 

525344669063 

53458836253 

38275293391 

264.758222487 

109  033459770 

200  399299969 

130  448227050 

564453719467 

16.387049764 

689  920420726 

124  595174314 

282  346963237 

302  451537677 

17  105846707 

298.005804703 

792704752916 

277.592370478 

488  070794225 

600  835486551 

77  330402620 

651  631129729 

186  384979015 

302942800237 

488  356688854 

282631614706 

12.094239822 

251  322489919 

639047414057 

315.553538531 

102  245119335 

121  065856465 

248  600567386 

120.292743494 

281  418670027 

6  512039507 

394  628880256 

329  862437764 

42  172683323 

275  823361976 

400377074953 

205  279868273 

89  822151642 

568897212277 

126  945901458 

413  319725633 

280829982419 

74621453179 

172  347938766 

138  081515648 

169  690890193 

295  648493372 

49  777202258 

271  146965769 

427  883569557 

120  358845277 

86  140070778 

106  453003164 

336696295308 

392  945334315 

439747178352 

145848786206 

357  670884257 

95  424425759 

158  394437081 

233.214847248 
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Row  Names 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 


V1 

V2 

V3 

78318115166 

233.292392917 

150.556841342 

400.322169035 

703.176710554 

131.985760999 

284.082094482 

700.398323845 

54.502213913 

673.813459742 

628.345709323 

27.810511050 

416.157442911 

309.802781935 

61  701311275 

502.727434228 

208.984463635 

7.000471997 

395.249167889 

181.615779362 

246.169600960 

532.055233005 

282.774380172 

51.530633142 

256.526079351 

265.073575146 

366.480188363 

791.013013348 

483.679778572 

383.122947523 

116.674678096 

203.429431816 

501.860086251 

216.019262031 

50.783576381 

346.268316228 

652.055312871 

297.710549841 

19.143684755 

57.143580125 

444.863422850 

53.260846841 

305.430041853 

64.522458661 

344.328881540 

134.875969833 

153.848996778 

74.573546508 

207.795329813 

197.474623424 

36.572283433 

55.079317773 

148.585617029 

153.485311966 

332.356937878 

183.167123218 

2.771444009 

33.292174799 

281.412899135 

521.055842681 

329.237256924 

624.659923834 

410.145804201 

134.999172015 

924.566774196 

959.801158180 

310.283780658 

196.949468983 

61.420855422 

145.927518733 

98.091275956 

105.358948583 

96.240971550 

448.577731504 

68.349212862 

504.457185330 

3.995392340 

54.819812986 

779.185987694 

63.103630524 

510.654006981 

64.295309451 

227.558484641 

681.717084473 

132.919226467 

80.090112060 

305.523448635 

225.006666940 

354.405022103 

115.791651829 

322.833702315 

15.741543588 

273.716914308 

552.256221658 

238.880493230 

103.236641646 

234.905588128 

267.774447132 

384.884639091 

241.821803683 

27.660947307 

35.559153707 

360.121368725 

54.400309193 

74.757012406 

109.375105806 

169.046509653 

306.549361255 

211.877188712 

54.427577933 

47.413167726 

328.701690509 

233.519929404 

411.204384691 

271.788277670 

132.777355320 

359.728088427 

54.715071977 

52.466467632 

253.216046834 

157.314280817 

49.880943701 

746.100150202 

93.608619218 

110.699690363 

573.390687287 

218.283430910 

434.282741182 

109.443120744 

1106.554817220 

178.424617471 

107.752314459 

179.703059410 

96.191152623 

63.123472845 

220  753490182 

358.476353401 

438.979104011 

97.987549536 

57.556717482 

152.972644520 

126.992796611 

173.350825362 

1004.747325381 

656.574251587 

250.152057881 

312.989740534 

189.330584986 

142.057479837 

63  629475350 

Column  1 

Column  2 

Column  3 

Column  1 
Column  2 
Column  3 

1 

0.012345873 

1 

-0.151985402 

0.027275574 

1 
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A         1          B          |          C 

0 

r 

H 

I 

K                                   M 

N                                                     Q 

1 

MAXIMUM  LIKELIHOOD 
ESTIMATOR 

Number 

of 

observati 

ons 

MODIFIED  MOMENT  METHOD 

2 

shape  parameter 

n 

. 

lailures 
limes  Ti 

shape  parameter 

3 

nitial  beta 

300 

mitral  beta 

4 

1  26241 

126241 

5 

final  beta] 

1  261028: 

fnal  beta 

6 

-1  9E-05 

1734948 

1288232 

7 

scale  parameter 

scale  parameter 

8 

landa 

landa 

9 

0  003418 

Ti 

LnTi 

TI«B4 

Ti"B6 

LnTi- 
Ti"B6 

Ti»Q6 

(Tl»Q6- 
mean)A2 

(Ln  Ti)"2 

0003394 

10 

theta 

4506347 

134  6203      763  602 

450  6347 

2.771444 

1  019368 

3  621402 

3  616305    3  686347 

3  717989 

2290781     1  039112 

theta 

11 

292  5569 

395  2643 

300  696'  201  4339    3952643 

3995392 

1  385142 

5746651 

5  735662    7  944705    5  955909 

2284012 

1  918618 

294  6595 

12 

467  906 

301  5104 

22  16672      467  906 

4  372291 

1  475287 

6  439285 

6  426172    9  480448 

6  689317     22817% 

11  17522     2268263 

12  26648'    2264978 

2  176472 

13 

200  0048 

261  4197 

234  0198    200  0048      6  51204     1873653 

1064743 

10  6199        19  898 

3  510574 

14 

62.18725 

351  5663    293  9545 

62  18725    7000472     1945978 
""  188  4206'   12  094241   2  492729 

11  66534 
232626 

32  44627 
34  13487 
36  03582 

1163401     22  63953 

3  786829 

15 

188  4206     530  4414     446  4993 

2318261 
32.32293 

57  78798    24  80917     2227382 

6  213699 

16 

517  0458:   212  6503    4  372291 

517  0458    15  74154:   2  756303 
*  i04  4066*  16  38705!   2  796491* 

89.09181   3483971 

2197542    7  597208 

17 

104  40661   91  44069 

286  231 

34  00322    95  089711   36  69092 

2192057 

7  820364 

in 

1855807    1896589 

4614523    18 

35  89471     1019202    38  77717     2 

8  062308 

19 

2788055 

142  6871      63  7859    278  8055,    19  14368'   2  951973 

41  53768 

41  36859 

122  1189    44  82767 
*  124  5595~*45  54801 

2168030 

8  714144 

20 

2423741 

140  0031:   57  57154!  242.3741     19  38205    2  964348 

42  19167 

42  0192 

2165909    8  787356 

21 

191  6384 

35  9825 1   135.6856 

191  6384 
r*1156514 

22  166721   3  098592 

49  98356 

49  77001      154  2169    54  14712      2140672 

9601272 

22 

115  6514     807  5153,   56  21773 

27  66095 

3  320022 

66.104 

65  80144 

218  4622    72  02099     2088689 
220  3101     72  52305      2087238 

1 1  02254 

23 

318  7385 

220  5115 
121  6994 

96  95887 

318  7385    27  81051 

3.325414 

66  55654 

66  25042 

1 1  05838 

24 

389  324 

3919678       389  324     33  29217     3  505322 

83  52575 

83  12216         291  37,   91  43863      2032940 
90  32207       322  558"  9953712*    2009912 
9168017"     328  493     1010663      2005578 

12  28728 

25 

108  5496 

833  6427,587  01011108  5496    3555915    3571198 

90  76888 
9213521 

94  04576 
99  60763 

12  75345 

26 

6472831 

270  0064!   600  5939    6472831 

35  9825     3  583033 
'36  57228    3  599291 

1283812 

27 

135  2779    235  2606 

294  6109    135  2779 

93  57919]   336  8187 

1032054      1999524 

12  95489 

28 

167  5656     19  38205 

80  68715    167  5656    38  27529    3  644805 

99  10723 

361  2265 
4190714 

109  4375 
123  9987 
126  204 

1981938'     132846 

29 

620  4171 

4394084     348  3392    620  4171     4217268     3741773 

112  5787 

111  9981 

1941151      14  00086 

in 

549  6525 

1472637 

402  3291 

5496525    42.75377 
160  8586'  43  94084 

3  755457 

114  5404 

113  94761  427.9253 

1935011>   14.10346 

31 

160  8586 

198  6978 
2164074 
4275377 

131  9361 
272  5643 

3  782844 

118  5697 
125  6921 
130  5185 
138  787 
139  1522 
142  3386 

117  9516    446  19251   130  7361 
125  0288!  478  7419    138  7547 
129  8244*  500  97941    144  1938 

1922423 
1900251 

14  30991 

32 

269  7486 

269  7486    46  01891     3  829052 

1466164 

33 

119  8829 

568  1059     119  8829 

4741317         38589 
49  7772     3  907557 

1885285    14  89111 

34 

290  3671 

46  01891     470  1354    290  3671 
5904016    4425726    1115262 

'  112  0684*  799  0738     3158697 

1380396^  5393977^  1535213     1859758 
138  4025*  541  1038*  153  9337*    1858633 

15269 

35 

1 1 1  5262 
315  8697 

49  88094     3  909639 
*  50  78358"  3  927573 

15  28528 

36 

141  5682 

4  556  0193 
*  568  4S93 

157  5314      1848837 

15  42583 

37 

258  4736 

740  0128,   730  3911 

258  4736 
936  0746 

5153063    3  942176 
'"52  46647,   3  960174 

144  987 
148  3189 

144  1994 

160  523 

1840710     15  54076 

38 

936  0746 

525  3447!  53  45884 

147  5095'  584  1633:    164  2883 

1830507     15  68298 

39 

38  27529 

264  75821109  0335    38  27529 

53  26085    3  975201 

151  1594 
151  8691 

150  3314 

151  0365 

597  5976     167  4997 

1821828     15  80223 

40 

200  3993 

130  4482    564  4537 

'  689  9204*  124  5952 

302  4515*  17  10585 

200  3993    53  45884     3  978912 

600  9608     168  3022 

1819662     1583174 

41 

16  38705 

16  38705)  54  40031 !     3  99637 

155  2534 
155  3516 

154  3984 
154  496 

617  033    1721302 

1809349 
1809050 

15  97097 

42 

282  347 

282  347 

54  42758     3  996871 

617  5005     172  2413 

15  97498 

43 

298  0058 

792  7048    277  5924    298  0058 

54  50221 

3  998241 

155  6206 

154  7632    618  7806     172  5457      1808232 

15  98593 

44 

488  0708 

600  8355       77  3304    488  0708    54  71507 

4  002139 

156  3882 

155  5258    622  4358 

155  901 3*"624  2369 

156  8325;  628  7062 

173  4143      1805896 

16  01712 

45 

6516311        186  385     302  9428 

6516311     54  81981     4  004052 
*  488  3567*  55  07932    4  008774 

156  7663 

157  7037 
161  8296 

i     173  842:    1804747    16  03243 

46 

4883567 

282  6316    12  09424 

174  9029     1801898 

16  07027 

47 

251  3225 

639  0474 

315  5535    251  3225 

56  21773     4  029232 

*  57  14358    4  045567 

*  57  55672'  4  052771 

160  9311     648  4288     179  5736      178938C 

1623471 

48 

1022451 1   1210659 

248  6006    102  2451 

165  2014 

166  7106 

164  2805    664  6077     183  3925 
"  165  7796*  671  8669*   185  1023 

1779178     16  36661 

49 

120  2927 

281  4187      6  51204    1202927 

1774619     16  42495 

50 

394  6289 

329  8624 
400  3771 

42  17268    394  628S 

57  57154!  4  053028 

166  7648 

1658335:  6721277    185  1637     1774456!   1642704 

51 

2. '5  8234 

205  2799 

275  8234     59  04016    4  078218 

172  1531 

1711857     6981326     1912708      1758223     1663186 

52 

89  82215 

568  8972     126  9459 

8982215    6142086    4.117749 

180  9624 
182  0061 

179.9356    7409299 

201  2637      173182; 

202  4484      172870; 

16  95586 

53 

4133197.        28083     74  621451413  3197    61  70131     4  122305 

180  9723    746  0232 

16  9934 

54 

1723479     1380815     169  6909     172  3479    62  18725       413015 

183  8176 

182  7715     754  8737     204  504" 

1 72330; 
1713103 

1705814 

55 

2956485       497772       271  147:   295  648sl  63  10363     4  144778 

187  2436 

186  1743     7716513     208  395 

17  17919 

57 

106  453     336  6963    392  9453       106  453]   63  62948 

4  153077 
4 155532 

189  2155 

188  1328*      78133*210  6348*    170724! 

17  24805 

I                  i 

58 

439  7472     145  8488     357  6709    439  7472*     63  7859 

189  803 

1887162    784  2163*2113021      170550; 

17  26845 
17  33462 

59 

95  42443    158  3944 

233  2146 

150  5568 

I  95  42443!  64  29531 1  4  163487 

191  7185 

1906187    7936385    2134785     169982; 

I                 I 

60 

1   78.31812]  2332924 

78  3181. 

64  52246    4  167013 

192  574 

1914683     797  851 

214  4506      169728E 

230  9735       165450' 

17  364 

61 

400  3222 

|  703  1767 

131  9858     400  322. 

68  34921 
|  74  5735S 

4  22463 

I     207.103 

205  8976      869  841 

178475 

62 

284  082 

700  3983:   54  50221     284  082 
628  3457     27  81051     673  813 

4  311786 

231   1907 
231  3782 

229  8173      990  923    258  418' 

230  0035*  991  8735'   258  632: 

1584659       18  5915 

63 

673  813! 

1  74.62145 

4  31242! 

|    1584120    1859704 

64 

416  1574    309  8028I  6170131     416  1574'  74  75701 

4  314242 

I  231  9089 

230  5305 

994  564" 

1   259  23771    1582597    18  61269 

70 


A 

B                 C                 0                 E                 F                 G                H                 1 

J                  K                   L                  M                  N                  O 

P 

Q 

R 

65 

'  5027274 

208  9845]  7  000472 
181  6158    246  1696 

502  7274]     77.3304 

4  348087 
4  360779 

242  032]  240  5822]   1046  0721  270  7903     1553664 

18  90586 

66 

3952492 

3952492    78  31812 

245  9411     244  46361    1066  052       275  254 

1542556 

1901639 

67 

532.0552 

262  7744 

51  53063    532  0552 

80  09011     4  383152 

252  9866       251  459 
255  3697    253  8251 

1102  183 

283  3029 

1522627    19  21202 

68 

256.5261 

265  0736 

366  4802 

256  5261 
791  013 

80  68715    4  390579 

1114  439 

2860265       1515913'    1927719 

69 

791  013 

483  67981    383  1229 
203  4294'  501  8601 

85  14007 

4  455975 
4  497832 
4  515691 

277  3467    275  6442 
292  3959:  2905843 
299  063.   297  2027 
308.0416    3061155 

1228264.   311  1667      1454639]    1985571 

70 

1166747 

1166747 

89  82215 

1306  999]  328  4058* 
1342  075    336  04*89" 

1413352]  20  23049 

71 

216  0193 

50  78358     346  2683 

216  0193 

91  44069 

1395238]  20  39146 

72 

652  0553    297  7105 

19  14368 
53  26085 
344  3289 
74  57355 
36  57228" 
153  4853 

652  0553    93  60862    4  539122 

57  14358    95  42443    4  558335 

30543T96  19115*  4566337 

1389  496!   346  3474!     1 371 01 5T  20  60363 

73 

57  14358    444  8634 

315604    3136224!    1429596    355  0263     1350766 

20  77841 

20  85144 

7-1 

305  43'  64  52246 

3188087'   3168034'    1446631 

358  7054     1342227 
358  9447J    1341673 

75 

134876      153849 

134  876    9624097    4  566855 
207  79531  96  958871  4  574287 

3190171     3170103       144774 
322  0242*319  9952*      1463  75 

20  85617 

76 

207  7953.   197  4746 

362  3977 

1333686 

20  9241 

77 

I  5507932,   1485856 

55  07932 

97  98755      4  58484 

98  09126!  4  585898 
1022451*1*4  627373 

326  3432 
326  7794 

3242823     1486  782 

367  3583 

1322253 

21  02076 

78 

!  332  3569    183  1671 !  2  771444 

332  3569 
*33  292ff 

324  715*2'    1489  111     367  8593 

1321101 

2)  03046 

79 

3329217    281  4129 

521  0558 

344.3447|       34215|   1583256 
348  5656^  346  339*3*    1605  984 
353  5599*  351  2963*   1632  928 

388  0482      1275098 

21  41258, 

81 
82" 



134  9992    924  5668 

959  801 2 
61  42086" 

134  9992 

104  4066 

4  648293 
4  657373 

398.6482I    1251272]  21  60663 

3102838    196  9495 

3102838 
1459275 
9624097 

105  3589 

357  636    355  34181    1654  959'   403  3386 
3623306*   360*0011*   1680378    408  7422* 

12408001   21  69112 

83 
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APPENDIX  E.  FUNCTION  "tocompare' 


funcuon(v,  betamle.  thetamle,  betammm,  thetammm) 

{  #    This  function  calculates  the  kolmogorov-smirnov  goodness  of  the  fit. 

#  and  plots  of  the  empical  distribution  and  theoretical  distribution 

#  for  a  set  of  data  and  two  sets  of  parameters  of  weibull  distribution. 
win.graphO 

a  <-  ks.gof(v,  dist  =  "weibull",  shape  =  betamle,  scale  =  thetamle) 

b  <-  ks.gof(v,  dist  =  "weibull",  shape  =  betammm,  scale  =  thetammm) 

par(mfrow  =  c(l,  2)) 

cdf.compare(v.  dist  =  "weibull",  shape  =  betamle,  scale  =  thetamle) 

cdf.compare(v,  dist  =  "weibull",  shape  =  betammm,  scale  =  thetammm) 

hst(a.  b) 

} 
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tocompare(a20. 1,0.973,251.71,0.9857,253.14) 
[[1]]: 

One-sample  Kolmogorov-Smirnov  Test,  hypothesized  distribution  =  weibull 

data:  v 

ks  =  0.0546,  p-value  =  0.994 

alternative  hypothesis:  True  cdf  does  not  equal  the  weibull  Distn.  for  at  least  one  sample  point. 


[[2]]: 


One-sample  Kolmogorov-Smirnov  Test;  hypothesized  distribution  =  weibull 

data:  v 

ks  =  0.0509,  p-value  =  0.9977 

Alternative  hypothesis:  True  cdf  does  not  equal  the  weibull  Distn.  for  at  least  one  sample  point. 
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Empirical  and  Hypothesized  weibull  CDFs       Empirical  and  Hypothesized  weibull  CDFs 


D 


-solid  line  h  Ihs  eni  pi'  >zal  d/ . 

MLE 


valid  line  is-  I  he  em  pi'kal  d/. 

MMM 


Figure  5.  Kolmogorov-Smirnov  Test;  Hypothesized  Distribution  =  Weibull.  (a)  ks  = 
0.0546,  p-value  =  0.994,  (b)  ks  =  0.0509,  p-value  =  0.9977 
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tocompare(al  00. 1,1.261,292.56,1.2882,294.66) 

One-sample  Kolmogorov-Smirnov  Test;  hypothesized  distribution  =  weibull 

data:  v 

ks  =  0.03 1 9,  p-value  =  0.9207 

Alternative  hypothesis:  True  cdf  does  not  equal  the  weibull  Distn.  for  at  least  one  sample  point. 

[[2]]: 

One-sample  Kolmogorov-Smirnov  Test;  hypothesized  distribution  =  weibull 

data:  v 

ks  =  0.0382,  p-value  =  0.7741 

Alternative  hypothesis:  True  cdf  does  not  equal  the  weibull  Distn.  for  at  least  one  sample  point. 


78 


Empirical  and  Hypothesized  weibull  CDFs       Empirical  and  Hypothesized  weibull  CDFs 


valid  line  ti  (.he  em  pi'kal d/. 
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solid  line  ivlheempirioldS. 

MMM 


Figure  6.  Kolmogorov-Smirnov  Test;  Hypothesized  Distribution  =  Weibull.  (a)  ks 
0.0319,  p-value  =  0.9207,  (b)  ks  =  0.0382,  p-value  =  0.7741 
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